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|. INTRODUCTION 


|. General Problem of Turbulent Transport 


In the transfer of momentum, mass, and heat from turbu- 
lent fluid streams to solid boundaries, the major resistance 
has been thought to reside in a laminar sublayer adjacent 
to the surfaces. The large number of semi-empirical and 
theoretical equations available for such cases involve the 
concept of molecular transfer through the laminar sublayer, 
with suitable corrections to account for the resistances in 
the transition layer and in the turbulent core. 

In another class of transport problems there are no 
solid boundaries, and the entire resistance to transport of 


momentum, mass, and heat is in the turbulent fluid. For 
instance, in a gas burner a combustible gas issues from a 
nozzle and entrains and mixes with air. The position of 


the flame 1s determined by the design of the combustion 
chamber, the turbulent transfer properties, and the rates 
of the chemical reactions involved. Similar considerations 
apply to the performance of combustion chambers using 
atomized or powdered fuel. In general the problem of 
transport in free turbulence is common to such operations 
as atomizing, spray drying, combustion, and jet injection. 

Combustion in a turbulent flame involves the turbulent 
transport of heat, mass, and momentum, the generation of 
heat, and the formation of the products of the combustion 
reaction. The local rates of heat release and mass conver- 
sion are functions of the temperature and of the composi- 
tion, which in turn are governed in part by the transport 
properties. Experiments in the complex situations just 
mentioned would probably not yield an understanding of 
the basic physical principles underlying the process. 
Instead it seems best to study simple cases and to synthe- 
size the general case from the results of experiments 
on its elements. 

Accordingly, the first case chosen for study was the 
transport of momentum in an isothermal, turbulent, free 
Heteot air discharging at subsonic velocities into air 
of the Same density. The study of this case involves 
several theoretical treatments and the presentation of 


much experimental data. 
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The next steps would be to introduce temperature gradi- 
ents, then density gradients, then both simultaneously, 
and then finally to include a chemical reaction involving 
heat release. 


2. Method of Investigation 


The purpose of the work reported in this bulletin was 
to obtain and correlate information on mixing in fluid 
streams. Particular emphasis was placed on the measurement 
of momentum transport in turbulent flow, with the object of 
correlating this information, together with other investi- 
gators’ measurements of mass and sensible heat transport, 
in a form which will have general engineering utility. 

Consideration was given first to turbulent mixing in 
free jets. Of the available theories of turbulent trans- 
port, that of Reichardt '58!*" was chosen as involving the 
fewest objectionable assumptions and as offering most 
promise of extension. It was discovered that the assump- 
tion embodied in Peichardt’s formula for turbulent shear 
stress also linearizes one of the equations of motion when 
applied to free jets. This important result makes possible 
the mathematical solution of problems involving complicated 
boundary conditions by the method of superposition of 
particular solutions of the linear equation. The pre- 
dictions resulting from the mathematical solutions are 
compared with the results of experiments performed in 
this laboratory. 

Reichardt’s hypothesis was also extended to cover the 
cases of the transport of mass and heat in free jets. The 
resulting equations are also linear, and their application 
to the correlation of experimental data is considered. 

The current theories of turbulent flow as applied to 
free jets are reviewed. To compare the assumptions and 
range of application of each of these and to clarify 
certain statements found in the literature, it seemed — 
desirable to develop the theories on a more comparable 
basis than was done by the original authors. In some cases 
this procedure required slight deviations and extensions 
from the treatment presented by the original authors. Some 
unstated assumptions which seem implicit in their original 
treatments are also mentioned. 


3. Experimental Program and Description of Systems Studied 


The experiments reported herein dealt with the transport 
of momentum in isothermal, turbulent, free jets of air 
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discharging at subsonic velocities into stagnant air 
The nozzles used were designed to provide a flow of face 
form velocity over the discharge opening and a low level 
“of turbulence. 

In such aeyet there are several well defined regions. as 
shown in Fig. 1. The nonturbulent core is a cone téhapad 
region having its base on the nozzle opening and its apex 


———— 
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Fig. |. Zones of Flow and Momentum Flux Profiles in Free Jets 


on the axis of flow about 4 nozzle diameters from the 
nozzle. In this region the intensity of the turbulence is 
low; the velocity is substantially uniform and equal to the 
Svelocity at the discharge. Surrounding this core 1s a 
wedge-shaped annular region with its apex at the lip of the 
nozzle. Here the flow is highly turbulent and the velocity 
decreases rapidly away from the nonturbulent core. 
In the region of flow more than 10 nozzle diameters from 
the nozzle the turbulence is said to be fully developed, 
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that is the profiles of the velocity versus the radius are 
similar at successive sections taken perpendicular to the 
axis of flow. The relative intensity of the turbulence may 
become 20 percent at some points in this region; 1t varies 
in a complicated manner.‘*°!’ Between the nonturbulent 
core and the region of fully developed turbulence is a 
transition region. 

The many types of mathematical functions that have been 
proposed for fitting the profiles of velocity in the region 
of fully developed turbulence include power series, |7°! 
trigonometric series, ‘°*) and the probability function. ‘?? 
This last is the most popular and involves the fewest 
arbitrary constants, ‘9°! 

The locus of points in the jet at which the momentum- 
flux density is one-half the value of that on the axis of 
flow defines the half-momentum-flux angle. This is an 
angle of approximately 10 deg, whose vertex lies near but 
not at the nozzle discharge. 

Two types of single free jets have been considered-- 
the two-dimensional jet which issues from a long, narrow 
slot, and the round jet which issues from a circular 
nozzle. The latter type was used in the experiments 
reported in this bulletin. The results of these experi- 
ments indicate that all such free jets are essentially 
similar when compared on a dimensionless basis involving 
ratios of velocities and ratios of dimensions in the jet 
field and the nozzle diameter. This finding conforms with 
general jet theory. 

Measurements were also made in the turbulent flow field 
adjacent to a nozzle discharge, in the field formed by the 
confluence of two parallel jets, and in the field formed by 
the discharge of a free jet parallel to a plane surface. 
Advantage was taken of the linearization of the equation of 
motion by Reichardt’s hypothesis to predict the distribu- 
tion of momentum in these several cases of free turbulent 
flow with complicated boundary conditions, The predictions 
are compared with the experimental results in Section 19. 

Although the experiments were restricted to the trans- 
port of momentum, the extension of Peichardt’s hypothesis 
to heat and mass transfer is applied to correlating heat 
and mass transfer data presented by other investigators. 

The experimental data were obtained by measurements with 
a total-head impact tube which was mounted in a special 
traversing mechanism permitting precise measurement of the 
position of the probe relative to the discharge nozzle. 
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hci ation 


Ay Discharge area of flow nozzle 

a Distance of equivalent point source from 
nozzle 

ay Proportiona lity constant in Eq. 43 

) ) 0) 6 


Jet breadth parameter 
Proportionality constant in Fq. 46a 


0 
b Term in Eq. 46a to account for transport 
1 : ; : 
of heat or mass in the x-direction by 
the longitudinal fluctuations of velocity 
Oy Term in Eq. 44a to account for transport 
of momentum in the x-direction by the 
longitudinal fluctuations of velocity 
ee ove once von 3 
es Spreading coefficient defined by Eq. 58 
Cy Impact tube flow coefficient defined by Eq. 76 
C, Specific heat of fluid 
Cy Tollmien’s spreading coefficient defined by 


Eq. 24 
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D Distance between two parallel nozzles, or 
twice the distance between a nozzle and 
a plane parallel to its axis 


D Mass Diffusivity 


Generalized symbol for molecular diffusion 
constant of heat, mass or momentum 


d Diameter of flow nozzle 

E Thermal expansion factor (see reference 3) 

|e nt Fraction of flux (mass or heat) transported 
by mean-flow components 

I Constant used in Eq. 66 

hy Pressure drop across standard flow nozzle, 
in. of water 

Likes a: ke Unit vectors in the x, y, z directions, or 
in the x, r, ® directions resfectively 

K Flow coefficient in Eq. 73, includes nozzle 
diameter 

Kp Strength of a free jet, Eq. 47 

K Term to correct for the nonuniformity of 
the velocity distribution at the nozzle 

Ky Term to account for deviation from a uni- 
form concentration or temperature dis- 
tribution at the nozzle 

k Thermal conductivity, heat units per unit 


area per unit temperature difference per 
unit length 


eee Ly: Lj, Vector mixing length for momentum, heat 
and mass, and vorticity 

Soph ad  r ae Scalar components of L' in cylindrical 

~ coordinates 

ea Lys ar Magnitude of L' under conditions of iso- 
tropic turbulence 

lp Pe, Prandtl’s mixing length defined by Eqs. 20 

lp Taylor’s mixing length defined by Eq. 28 

M Molecular weight 

N Unit vector directed outward normal to 


surface 
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Pressure 


Stagnation pressure 


Pressure upstream from standard flow nozzle, 
psia 


Pressure recorded by impact tube manometer 


Generalized symbol for local rate of gen- 
eration of heat, mass, or momentum 


Local rate of generation of sensible heat, 
heat units per unit of volume per unit of 
time 

Local rate of generation of i-th component 
of fluid, mass units per unit volume per 
unit time 


Discharge capacity of standard flow nozzle 
in cu ft per sec measured at 60 deg F and 
14.7 psia 

Universal gas constant. Also a vector dis- 
tance 


Correlation coefficients defined by Eqs. 18 
and 37 

Polar coordinates in any plane perpendicular 
to axis of flow 


Unit vector in any given direction 

Polar coordinates in the plane of the nozzle 
erifice 

Absolute temperature, deg R 

Temperature upstream from standard flow nozzle, 
deg P 

Difference in temperature between that of a 
point in a free jet and the ambient tem- 


perature, t = (T <ie+) 

In both rectangular and cylindrical coor- 
dinates: Instantaneous velocity paral- 
lel to x-axis 

Velocity vector (components u, v, and w) 

In rectangular coordinates: Instantaneous 
velocity parallel to the y-axis. In 


cylindrical coordinates: Instantaneous 
velocity in the radial direction 
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w In rectangular coordinates: Instantaneous 
velocity parallel to the z-axis. In 
cylindrical coordinates: Instantaneous 
angular velocity 

vw, Concentration of i-th component in fluid, 
mass units per unit mass of fluid 

ARON pas Cylindrical coordinates 

<a Siete Rectangular coordinates 


Any scalar or vector function of the space 
and time coordinates 


Oo, Bor By Constants in Eqs. 45 and 46 


Op, Op, Op "Eddy diffusivity" of heat or mass defined by 
equations following Eq. 27 

SBeeP tT "Eddy viscosity" defined by Eqs. 15, 35 
and 43 

€ Base of the natural logarithms 

a z or ¢ component of w, the vorticity vector 

n y or r component of w 

6 Time 

A Proportionality function defined by Eq. 50 

bh Viscosity of fluid, force times time per 

unit length squared 

& x component of w 

p: Density of fluid, mass per unit volume 

Py Density of fluid upstream from standard flow 
nozzle, lb per cu ft 

2 Summation operator, sigma 

o Surface of an arbitrary volume, 7 

a Arbitrary volume in a flow field 

$ Generalized symbol for molecular transport 
potential. Also compressibility factor 
in Eq. 73 (see reference 3) 

y Generalized symbol for concentration of 
sensible heat, mass, or momentum 

w Vorticity vector (Components €,7,Z) 


V Vector differential operator, del 
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Line above symbol indicates a time average 


ma IGE aNG 
u2 = Hal 1 u2d@ 
Oy 
1 64Aé 
Sea taal af ud@] 2 
1 


primed symbol denotes a fluctuating component 


The following subscript notation is also used: 


a Ambient conditions 
dA Conditions at an elemental point source jet 
h Sensible heat 
j-th point source jet 
m Momentum 
s Streamline of flow 
t Temperature 
w Mass 
Seni Conditions at points (x,r) in jet, any ¢ 
x0) Conditions at point x on axis of flow 
0 or 0,0 Conditions at discharge of nozzle 
vA Radius at which a flow variable has one-half 


the value on axis of flow 


11. MATHEMATICAL BACKGROUND AND PREVIOUS 
CONTRIBUTIONS 


6. Transport Equations 


A generalized differential equation for the flux of mass, 
heat, or momentum may be obtained by equating the rate 
of accumulation in a volume 7 to the rate of generation in 
the volume plus the rate of transport through the surface 
of the volume co. The rate of transport may be considered 
as that resulting from molecular motion plus that from 
the bulk motion of the fluid. Let % denote the concentra- 
tion of mass, heat, or momentum, as the case may be; then 
the rate of accumulation within the volume is 


2_f bdr = f Qdr + f DVO-Ndo - f{ BV" Ndo (1) 
Mela T a o 


where the first term on the right represents the rate of 
generation within the volume r and the second and third 
terms represent the net rate of transport across the bound- 
aries of the volume by molecular processes and by motion of 
the fluid, respectively. For the three cases the general- 
ized variables, ~, Q, D and ¢, are interpreted as shown in 
the table below. 


Interpretation of Generalized Variables 


p Q D p 
Mass flux pw; 0) F [ye pu ; 
Heat flux pC, 0, k v 
Momentum flux pVv's -VP’S LL Ves 


Provided a volume exists which on the one hand is 
suffic1ently small that all the flow parameters differ from 
their space-mean values only by second-order differences 
throughout the volume and which on the other hand is 
sufficiently large that the statistical properties of 
velocity, pressure, etc., are defined, Eq. 1 may be trans- 
formed by the divergence theorem, to 


ae 


56° Q+V'DVE - V'PV (2) 


16 
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For application to free turbulent jets it is now assumed 
that the pressure is uniform throughout the flow field, 
that there are no chemical transformations, no heat release 
or absorption by a chemical or physical process, and that 
the amount of sensible heat formed by conversion from 
kinetic energy is negligible. Under these conditions the 
generation term Q in Eq. 2 vanishes. [Furthermore, under 
most conditions of free turbulence, the molecular transport 
term is negligible. Hence, for the quasi-steady state of 
steady turbulent flow (i.e., no change with time in the 
mean flow parameters) 


EN Verws 0 (3) 


The bar denotes the operation of time-averaging: 


Z fs Lim 1 (at? 
60 6 ~O, 


Zd0 


where Z is any fluctuating variable and where it is sup- 
posed that the value obtained for Z is independent of the 


arbitrary reference time 9). Hence 
8. 
= 0 


A special case of Eq. 3 is now developed. Substituting 
pw, for ¥ and summing over all components present, 


V-(Zow,)V = 0 
i 


But by definition 


whence 
V- pV = 0 (4) 


which is the familiar equation of continuity. Introducing 

the concept of resolving the flow parameters into time-mean 
(30) 

and fluctuating components, 3 ret 


Vos ee oes Veto V LAvector. addition) 


= | 


u= 5S Uh 
Fes, 9 y) 

vo = VU +t -U 
me ’ 

w=w +w 


AS) 
W 
1 
: 
d 
. 
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where, by definition 
VE c= 


u' = 0, etc. (6) 


tt 


For the special case of an incompressible, homogeneous 
fluid in isothermal flow (flow with uniform density), T and 
e do not fluctuate with time. 


p' = 0 
VipV = 0 (7) 
Expanding Eq. 4 for the instantaneous case (i.e., Eq. 4 
with bar omitted) 
VipV + VepV' + Vea'V + Vep'V' = 0 (8) 


Applying Eqs. 6 and 7 for the case of flow with uniform 
density, 


pv:V" = 0 (9) 


For use in later arguments, consider now the product 
ZV:pV, where Z is any fluctuating property of the fluid. 
Let Z = Z + Z' and expand the above product by the aid of 
this relation and Eq. 5. Now applying the same arguments 


which led to Eq. 9, it follows that 
ZV:pV = 0 (10) 


for incompressible flow with uniform density. 


7. Phenomenological Theories 


a. Prandtl’s Model ‘5? 


In developing the mathematical analysis of isothermal 
free jets, the following restrictive assumptions are 
frequently made. 


1. The mean pressure Pand the mean density p are 
uniform throughout the flow field; hence their 
derivatives are zero 

2. The jet and ambient fluids are incompressible and 
have the same density and temperature throughout; 
hence po’ is always zero 


The second is an important assumption, for if p'’ assumes 
finite values, Eqs. 7, 9, and 10 are only approximate. For 
example if po’ is not zero, then the ideal gas law is not 
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satisfied by the mean components of pressure, temperature 
and density. The law may be written 


p = Ror 
Mu 


‘Substituting the mean and fluctuating components and 
obtaining a time average of the resulting expression gives 


—— R ane oan, 
Peet oT aie 
Pe tek ) 


Fluctuations in density and temperature due to pressure 
fluctuations are always of the same sign; hence the product 
pT’ never assumes negative values, and the mean value 
therefore cannot be zero unless p' is always zero. 

Consider a free turbulent jet issuing from a point 
source into an ambient fluid having the same density and 
temperature as the jet fluid, and choose a system of 
cylindrical coordinates (x,r,¢) with the jet discharging 
along the x-axis. For the x-directed momentum flux, Eq. 3 
becomes 


VeouVv = 0 Gis) 


Applying the two assumptions listed above, together with 
Eq. 8, this becomes 
V-Vu =- Veu'V' (125 

By introducing the concept of mixing length, Prandtl ‘9? 
deduced a relationship between the u’V’ product and the 
mean flow parameters and their derivatives. Though Prandtl 
originally treated the case where the turbulent motion 
as well as the mean motion is confined to two dimensions, 
his method is readily extended to cover motion in three 
dimensions, as follows: It is assumed that momentum 
is conserved, so that a parcel of fluid at the point 
(x ,r ,¢ ) having the mean momentum associated with the 
flow at that point is conceived to move without alteration 
of its original momentum along a vector path L’ to a mixing 
Point, (x,7,¢). Because of the turbulent motion, L’ is 
seldom parallel to the mean velocity vector at (x ,r ,¢ ie 
Further, where there are velocity gradients, the mean 
velocities at the two points will, in general, be different. 
The parcel of fluid, therefore, has an excess or deficiency 
of momentum, and upon mixing with its surroundings at 
(x,r,@) causes a fluctuation in the velocity at that point 


20 ILLINOIS ENGINEERING EXPERIMENT STATION 


proportional to the difference between the mean momentums 
at the two points. That is, for x-directed momentum 


pu' = pdu = plu, r,o- u, ,r ue 


The model may be analyzed in the following way. The 
increment, Au, may be expanded in terms of the derivatives 
of u in a Taylor series, and if all terms of higher order 
are neglected, 


Au. = u’ =- Ei Va (13) 


where Ee is the vector path along which mixing occurs to 
produce u’ The vector mixing length L’ is represented as 
a fae chat tan since it has significance only in connection 
with a fluctuating quantity such as u’. Substituting into 


Eq. 12 


V-Vu = Vel Lio va k (14) 
L' and V' may be expanded in terms of their scalar compo- 
nents. To obtain Prandtl’s result, it is assumed that the 
turbulence intensity and mixing length are isotropic, 1. 
which case the parallel products are all equal. 


Hed u! 2 Pe v! = it Ul = € (LS) 


where €, is defined by this relation; it is commonly called 
the eddy viscosity. 

It 1s reasonable to assume that to a first degree of 
approximation all cross products are zero. ‘?°) 


5 
- 
s 
J 


erie = Lo stoi" = A eae = Lew = Ly yu' = Lye’ = 0 
(16) 
Substituting into Eq. 14, t 
hae ES si 
V-Vu = Vee Vu (17) ¢ 


A relation between €, and the flow parameters is required — 
to solve Eq. 17. Let a correlation coefficient be defined 
such that 


ie ‘ 
R a ruv 
es (1! 2 yl2\A (18) 
ru ' 
In isotropic turbulence u'? = v'2 = y'2, Combining this 


and Eqs. 13; 15. and'18 


epee Rp Lieut vy ele (19) 


u 
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Again, in isotropic turbulence, it is assumed that 


eerie a Tp 3. 4 


xu’oru Li pat ru’ gu 


and 


Where 1’ is defined by this last relation. Hence 


Ep =e ip ead aire pale 


u 


or 


Ep = Lp *|Va (20) 


where 1p is the Prandtl mixing length and contains (impli- 


- u . - - . . 
citly according to Prandtl’s original derivation) the 


correlation coefficient Rp pubstitutine into Eq.) 17 
u 


V-Vu = V> p|val Vu ) 
(21) 
Imposing the condition of isotropy does not seriously 
_ limit the application of Prandtl’s result. For turbulent 
mehearing stresses to exist’it is not necessary that the 
mean-square fluctuating components differ from each other; 
“it is sufficient that the correlation coefficient in Eq. 18 
be greater than zero. A reasonable first approximation is 
that the difference between the mean fluctuating components 
is of secondary importance compared to the degree of 
correlation in determining the shearing stress. 

Equation 2] is a generalization of Prandtl’s early 
mixing length theory. His result is equivalent to the 
hypothesis that 


ey lp 2 |u| (22) 


which applies strictly when the turbulent motion 1s confined 
to two dimensions. The quantity \Vuloin Eqs... 20 cand 
21 reduces to du/dor if the terms du/dx and (1/r) (du/2¢) can 
be neglected. For axially symmetrical flow in free jets, 
where du/dd is zero and du/dx can usually be neglected 
compared to du/dor, Eq. 21 becomes 


Taw thee 7G) 2| du | (du (23) 
SW Yu (ec 
us r ao ltl p is | 
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which is the form given by Prandtl and applied by Tollmien 
to flow in free jets.‘’°’ The utility of Eq. 23 depends 


on the existence of simple relations between Lp and the 


boundary conditions for a particular case. For flow near a 
plane surface Prandtl supposed that lp) 1s proportional to 


u ix: 
the distance from the surface,‘?*’ and for flow in jets or 


wakes he assumed that 1p) was independent of r and was 


proportional to the breadth of the turbulent zone. ‘53! 

The spread of free jets is approximately linear with the 
distance from the nozzle discharge, ‘’°’ and the radial 
distribution of the momentum at successive sections through 
the jet perpendicular to the x-axis are similar. Hence, 


Tollmien ‘7°’ assumed 
lp = Cyx (24) 
u 
and fe 
r 
——tsrl= f(—) (25) 
u x 
x,0 
Corrsin ‘+?’ has shown that these two assumptions are not 
p 


independent, if it is assumed that the distributions of 
both the mean and fluctuating motions are dynamically 
similar at successive sections. = 

Tollmien'’°’ further assumed that the quantity du'?/dx 
is negligible compared to (l/r)(du‘v'/dor) and that the 
quantity ou/dx is negligible compared to (l/r)(0u/dr). 
These assumptions are equivalent to each other in view of 
the relation between u’ and Vu in Eq. 13. 

The assumption of similarity enabled Tollmien to define 
a stream function and then to separate the variables. He 
obtained a series solution which he presented in the form 
of tables of dimensionless parameters. From experiments 
performed at Gottingen, he estimated that C, is 0.156 
for axially symmetrical flow. When this value is used, 
the agreement between predicted and experimental mean 
velocities is fairly good. 

The principal defect in this treatment is that Eq. 24 
is incompatible with the fact that Rp must vanish on the 


u ; 
axis of a free jet because of the local isotropy of the 
turbulence. Hence, for lp to remain constant across a 


. - u . - . 
section, 1’, must become infinite at the axis -- an impos- 
sibility in the model assumed. Furthermore, 1’, cannot 
even become very large as compared to the dimensions 
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of the jet section, or the assumptions embodied in Bq 
therefore become invalid. 

Profiles are similar only from a point source jet, and 
the solution based on similarity woutd be expected to hold 
only in regions where the distance from the discharge is 
large compared to the diameter of the nozzle. Actually, 
Tollmien’s solution does not apply in the region within ten 
nozzle diameters from the nozzle. 

On the other hand, the application of similarity to the 
problem of two-dimensional turbulent flow around a corner 
leads to a solution that is valid in the immediate neigh- 
borhood of the corner. Tollmien solved this case also, and 
the solution was applied by Kuethe‘*°’ to the annular 
turbulent mixing zone adjacent to a circular nozzle. His 
method gave good agreement with measurements of velocity in 
this zone close to the nozzle but became impractical in the 
region where the turbulent mixing zones from opposite sides 
of the nozzle intermingle. 

The foregoing derivation of Prandtl’s theory may be 
extended to the general case of heat and mass transfer by 
fevaining Y in Eq. 3 so that instead of Eq. 12 there is 
obtained 


V-V¥ =- VW! (26) 
and corresponding to Eq. 13 there is 


yp! =- L'y*Vu 
From relations corresponding to Fq. 15 
oe = | ry” 


where $, 1s the eddy diffusivity of heat or mass as the 
case may be and corresponding to Eq. 17 1s 


V-Vy = V-8V¥ (27) 


By defining a correlation coefficient similar to that in 
Eq. 18 and proceeding as before, there results 
% fh Wit Piva \ oleh 
eee ae ply (vu)? | 
= i es 


The simplest assumption is that Rp. = Rey and Lp ‘ 


eee Then Eq. 27 becomes 


V-V¥ = Yelp? [|Vu] (Ve) ] 
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and by analogy to Fq. 23 
oop eb 2 [pp 2 out By (28) 
WOME sacle Me teva 


Jowarth'3°’ states that since Eq. 28 is of the_same 
form as Eq. 23 and since the boundary conditions of ¥ are 
the same as thcse of u, then the distribution: of ¥ must be 
similar to that of u, namely 


agrieeae re 
Tee Ata te (29) 


> 18 a solution of Eq. 23. 


where eee 1 ts 

Equation 29 is one solution of Eq. 28, since substitu- 
tion in Eq. 28 reduces that equation to Eq. 23. However, 
Hinze!34) measured the distribution of heat and mass in a 
free jet and found that they are not the same as the 
distribution of momentum; hence Howarth’s assumption about 


the equivalence of Rp and lp in the two cases is not valid. 
Hinze also investigated the assumption that 1, is dif- 
ferent from 1, , but was unable to obtain agreement between 


experimental and predicted temperature distributions. 

In genera!, both the theoretical and the experimental 
objections to the assumptions required for a solution for 
the free jet problem using the Prandtl mixing length are 
serious. In addition, it 1s shown later that there is 
reason to believe that the mean flow parameters (u, J, 
etc.) can be correlated only approximately on the basis of 
measurements performed with total-head impact tubes and 
thermocouples in highly turbulent streams. 


6b. Taylor’s Model 


Prandtl’s theory was based on the assumption that the 
momentum 1s conserved during the mixing process. Taylor 
investigated ' the corresponding assumption that the 
vorticity is conserved. In order to apply this concept, 
the vorticity 1s introduced explicitly into the equations 
of motion. Equation 1] may be expanded to 


pV “Va - +E ov = 0 (30) 


However, the second term is zero by Eq. 10. Partially ex- 
panding the first term into mean and fluctuating components, 
and imposing the condition of incompressible flow, 
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V-Vu =— Prous a ee os ou! 
Ox or =o (319 
Let a vorticity vector be defined as 
Q= Ei + wy + Ck 
where AS) 
ene 1 d(rw) 1 ov 
ror Ps op 
‘ Ik Gz Ow 
f ip Cla Ox 
Sy eA 5p (33) 
Ox Or 


and t, j, and k are unit vectors in the x, r, and ¢ direc- 
tions, respectively. From 7 and € the following equali- 
ties may be written 


Hg! ov’ ou’ ’ 
v St = 
Ox Or 
: ; ious ow! 
w'(-7' +—— - —) = 
r op Ox 


Adding these two equations to Eq. 31 and combining terms 
yields 
meet it se 1 ov’? Sail ct epee Ewe 
V-Vu =- — int! ee = Tie 
2.8% (34) 


At this point the assumption regarding the conservative 
transport of the vorticity is introduced. Taylor’s own re- 
marks ‘°°? on this subject are perhaps the most cogent: 


Unfortunately, the assumption that the turbulent com- 
ponents (of vorticity) are statistically isotropic does not 
give rise to any simplification in the analysis. No fur- 
ther progress can be made without introducing further 
simplifying assumptions. One such assumption is that the 
components of vorticity are transported unchanged by the 
turbulence just as the momentum is conceived to be trans- 
ported in the Momentum Transport Theory. A part of the 
fluid is conceived to leave a certain position with the 
yorticity components of the mean motion and to retain those 
components till it mixes with its surroundings after 
traversing the mixture length. This is clearly a very 
drastic assumption, because it is certainly untrue in 
detail, though it may be true when considered in relation 
to the effects on the mean motion. 
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With the above assumption, equations analogous to Eqs ea 
may be written 


&! rth Deve 
Tape een 
Csr nee rt Vo 


In isotropic turbulence, the cross products are presumed 
to be zero 


Oe ie ie ah = ee = ee = ul oe = ul Pere 0 
while the parallel products are assumed equal 
all" Ve = vilt iy ow ET as = €¢ (35) 


where €, as defined by Eq. 35 1s an eddy viscosity and in 


particular is the eddy viscosity associated with the vorti- 
city transport theory. Substituting these in Eq iioas 
and neglecting the term 0V'2/dx as is done in Prandtl’s 
momentum transport theory, 


The terms in parentheses may be expanded by Eqs. 32 and 33 
and the resulting terms regrouped and combined to give 

ae 3 
“V = - em Vv 
V ep Vu peeved 


x 


But by Eq. 7 V-V is zero for incompressible flow, hence 


which may be compared to Eq. 17. 
To evaluate the eddy viscosity Eq: Howarth'3°) simply 


took it to be of the same form as the €, of the momen- 
tum theory defined by Eq. 15. This result follows only 
if important assumptions are made. Let a correlation 
coefficient R, be defined as 

wits 

ERTS 15) as a 

Ry = eg (37) 
ry 
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If the Mea lelce Ueetaacuropieyom.. srw. Tniorder £6 
eliminate u'* it is necessary to use the momentum transport 
Beeory as embodied in Eq. 13. These substitutions in 
Eq. 37 lead to 


! 2 1 i y 
rll rn (L el 


! =e) aD wena 6 oD, 
oe ge rt ) 


where 1’, is defined by this equation. Hence 


nah et 2) un] 4 


or 


eer bae Via (38) 


where 1, is the Taylor mixing length and contains R,- From 


Eqs. 20 and 38 ep and €, are of the same form, provided 


that the momentum transfer theory may be used as well as 
the vorticity transport theory. From Eqs. 36 and 38 


V-Vu =-1,?|Vu| (V?u) (39) 


Because of the rotational symmetry, derivatives with 
respect to $ are zero; and if, as in the momentum theory, 
derivatives with respect to x may be neglected, then 


ae, Se ae nn (39a) 


This may be compared to Eq. 23. These procedures origina- 
ted by Taylor and carried forward by Howarth, making use of 
the vorticity transport theory, result in another equation 
of motion (Eq. 365) which differs from that obtained from 
the momentum transport theory (Eq. 17). The eddy viscosity 
which appears in both equations is the same function of the 
velocity gradient (Eqs. 20 and 38). 

Taylor’s theory suffers from the defect that it requires 
one to make all the assumptions that were made in deriving 
Prandtl’s theory (e.g., conservation of momentum during 
mixing) plus the assumptions regarding the VOrtLGLLy wat 
is uncertain whether some of these assumptions are not 


mutually exclusive. 
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If Eq. 21 is expanded in rectangular coordinates, then 
for the case of two-dimensional flow in the x-y plane, 
partials with respect to z are zero and, to the same degree 
of approximation as before, partials with respect to x 
may be neglected. Further, let us suppose that lp is 
independent of y. Then it follows that 


er ipo eee ne 
V-Vu =-21,2 —— ( 40) 


Proceeding similarly from Eq. 39, 


Be sa > du d7u 
VV lace 
oy dy? 


which is identical with Eq. 40 provided that lp = 2 lp. 


Abramovitch'?? has taken this result to be a general rela- 
tion between the two mixing lengths, but it holds only 
incidentally for two-dimensional flow, probably because in 
two-dimensional flow vorticity is conserved, ‘3°’ 

No such simple relationship exists for the axially 
symmetrical flow. 

While the momentum and vorticity theories give the same 
distribution for momentum flux in two-dimensional jets, 
they do not predict the same distribution in axially sym- 
metrical jets, as Howarth has shown. ‘3°’ He found that the 
momentum theory provides the best correlation, and he 
presented his solution in tabular form. 

The application of Taylor’s theory to the transport of 
heat and mass is not straightforward. Equations analogous 
to Eq. 32 cannot be written, because temperature and con- 
centration are not vectors. Although an equation analogous 
to Eq. 32 cannot be derived by these methods, one may 
simply assume, as Howarth did, ‘3°’ that 


Veve =< l 2 — — —( 


— (41) 
T Orr or ea: 


where it was also tacitly assumed that the mixing lengths 
for heat and mass are the same as for vorticity. A solu- 
tion of Eq. 41, obtained from the solution of Eq. 39a, 
provided a better correlation of data from axially symme - 
trical jets than did the corresponding solution of Eq, 23, 
However, in view of the complete lack of a phenomenological 
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Dasis for Eq. 41, this result can hardly be interpreted as 
a confirmation of the vorticity transport theory. 

The foregoing discussion of Prandtl’s and of Taylor’s 
theories states clearly the comparable basis on which each 
pests. Itmust be concluded that, the usefulness of each is 
limited to special cases. Neither theory offers promise of 
extension for use in the general treatment of mass, heat 
and momentum transfer in turbulent streams. 


8. Empirical Theories 


a. Boussinesq’s Hypothesis 


By analogy to viscous flow, Boussinesq'?°’ proposed 
that the u'v'’ product in the expansion of Eq. 12 can be 
expressed in terms of an eddy viscosity and the transverse 
gradient of the mean velocity, thus 


he oe oi eg 
V-Vu =- Ve(u eet “B54 oe (Uh W) k) (42) 
r 


For flow in a pipe, the term €, was supposed to be a func- 
tion only of the Reynolds Number, but was found to be a 
complicated function of the position in the pipe. For this 
reason Prandtl and Taylor abandoned the analogy to viscous 
flow and developed their theories along different lines. 
There is a formal resemblance between Eqs. 17 and 42. If 
the terms du/ox and ou'2/dx are neglected, the equations 
become identical. However, they are not based on the same 
fundamental postulates. 

Although Boussinesq’s hypothesis did not prove satisfac- 
tory for flow in a pipe, Hinze'?3’ felt that an investiga- 
tion of its applicability to flow in free jets would be 
worthwhile. He assumed €, to be a function of the velocity 
along the axis of flow and the distance from an equivalent 
point-source jet, and to be constant across a section of 


the jet. . Thus 
ein = dgu, o(xta) (43) 


where x is the distance from the nozzle and a is the 
distance from the nozzle to the equivalent point source. 
With these assumptions it is possible to solve Hinze’s 
equations to obtain 


ae VK, /8 d : 


_ a (44a) 


(x+a) (ag+y)” [1+ 


Ba, (x+a)? 
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where 6, accounts for the transport of momentum in the x- 


direction by the longitudinal fluctuations, u’?, <The 
quantity was evaluated from Corrsin’s data‘'*®’ and found 
to be approximately constant with an average value of 
0.00045. Hinze then determined the constants a and a, from 
his experimental measurements, obtaining 0.6 d and 0.00196 
respectively, where d is the nozzle diameter. K, 1s a term 
to correct for the nonuniformity of the velocity distribu- 
tion across the outlet of the nozzle. By inserting the 
numerical values of the constants and by making certain 
other approximations and assumptions, Hinze obtained the 
equivalent of 


Cee Sac 1 
— = ——___ {__———_} ( 44b) 
(40260) Pieee et ee 
0.0157(x+0.6a)2 


From his measurements, Hinze showed that the correlation 
obtained on the basis of the assumptions leading to Eq. 42 
is as good as that obtained by Tollmien and Howarth on the 
basis of assumptions regarding the Prandtl mixing length. 

The hypothesis of Boussinesq is readily extended to the 
transport of heat and mass. When this is done Eq. 26 
becomes 


oy 
j + Pw" k) 


VeVP =- Ve(piu’ i + 8, ao 
2 


Hinze first supposed 6, to be proportional to €p,, but found 
experimentally that this is not so. The ratio 5,/ep 


decreased with increasing distance from the axis of flow, 
and was not even approximately constant in the central 
portion of the jet. 

He then proposed the relation 


a tar ee (45) 


B 
: age ee 
a ries VK Y/ 8a f Bo 370/80 
POMUMUBys\l bce. ie By 5 22 
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eeewhich 6). accounts for the transport in the x-direction 


due to the fluctuating motion and Ky accounts for any 
deviation from a uniform distribution at the nozzle outlet. 
After evaluating certain of the constants and neglecting 
certain others, he obtained the equivalent of 


Gop 227 3 1.196 eas 


} 


2 
y (708d 0196+ 14-1? 
aaa Oe ShGr0 8d)” (46b) 


Hinze computed Bo/a5 bo betl.22 "and By/a2, to be 0.24. 


Equation 46 correlates very well Hinze’s data for the 
central part of the jet, but gives high values for the 
concentration and temperature in the region where the 
Patio of r/x is greater than 0.12. 

The principal disadvantage in the use of Hinze’s devel- 
opment stems from the complexity of the expressions and the 
large number of empirical constants introduced. Equation 
43 is certainly not true, and the success of Hinze’s corre- 
lation, as well as those of Tollmien and Howarth, seems 
largely due to the fact that any reasonable assumption 
regarding the dependence of u’v’ or }'V' on the mean flow 
parameters results in a solution that fits the data reason- 
ably well. As Corrsin has pointed out, ‘7°’ it does not 
seem possible to differentiate among these theories on the 
basis of mean-flow measurements made with a total-head 
impact tube. 


b. Reichardt’s Hypothesis 


An empirical approach to the problem of free turbulence 
was proposed by H. Reichardt.‘9°’ Noticing that experi- 
mental measurements in a free turbulent jet can be well 
correlated by probability functions of the type 


ary Kr =(r/ we (47) 


he investigated the conditions under which such functions 
become solutions of the equation of motion. Eqta tron 
as Eq. 47 have been used by several investigators. ©’ ~’ 
Kp, 1s a constant, depending upon the strength of the 
Hiet;) 6 1s a.parameter, depending only pn. *pUahd 9b Ls 
independent of r. 
Now for axially symmetrical flow and where density 
fluctuations may be neglected, Eq. 11 becomes 
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2 Oe ae (48 ) 
a sesttr uu) = 0 


[Teichardt showed that Eq. 47 is a solution of Eq. 48 
provided that == 


ah bu2 
uv =- Aso (49) 
where 
b db 
A Oo idy, (50) 


and 6-1's)a function of** ‘only. 

No phenomenological basis for Eq. 49 was proposed, but 
Reichardt remarked that the hypothesis must stand or fall 
depending on whether it permits a mathematical analysis of 
free turbulence and not on whether a mechanical interpre- 
tation can be found. Actually such a physical basis can be 
developed if one assumes that 

1. The mean flow is essentially one dimensional (x- 

direction) 

2. The intensity of the turbulence is small 

3. The radial velocity fluctuations are proportional to 

the mean velocity 
This 1s no more rigorous a theoretical basis than are the 
mixing length theories and hence is not presented_ here. 

In some of the preceding methods, the term du’?2/dx was 
neglected. This ommision may introduce considerable error. 
Peichardt’s procedure avoids the necessity of neglecting 
these longitudinal velocity fluctuations, since these are 


2 (20) 


contained in u*, Corrsin reports experimental values 


Of the ratio (u' 2/u2)% at various points in a free jet. 
Along the axis of flow it increased from a negligible value 
near the nozzle to somewhat more than 0.20 at 12 nozzle 
diameters from the nozzle, and remained substantially 
constant thereafter. This indicates that it is probably 


inadmissable to neglect u'2 and its derivatives. 


The Reichardt hypothesis merits further investigation. 
Equation 47 provides remarkably good correlation of ex- 
perimental measurements and contains only one arbitrary 


constant. Equation 49, because of its simplicity, pro- 
mises to be useful for the analysis of more complicated 
cases. An extension to cover these cases is described 


in the following chapter. 


lll. THEORETICAL ANALYSIS OF FLOW IN FREE JETS 


9. Generalization of Reichardt's Hypothesis 


The methods of Prandtl and Taylor neglect the contribu- 
tion of the fluctuating motion to the transport of momentum 
parallel to the axis of flow in a free jet. This imposes 
an unreal conservative relation on the square of the mean 
velocity u. Peichardt’s treatment uses the momentum flux, 
which 1s a truly conservative quantity and makes no assump- 
tion regarding the distribution of this between the mean 
and fluctuating components. The resulting equation contains 
only one arbitrary constant, and several of its investiga- 
tors'?*58)°9) have shown that it correlates the data well. 

Reichardt’s procedure can readily be generalized to in- 
clude the transport of heat and mass as well as momentum in 
free jets. This is especially useful because Reichardt’s 
hypothesis (Eq. 49) linearizes the equation of motion for 
x-directed momentum in a free jet and offers promise of 
making possible the mathematical solution of cases of free 
turbulence involving complicated boundary conditions. 

For axially symmetrical flow Eq. 3 takes the form in 
cylindrical coordinates, 


ier Lee Se 
= Soper, = 0 C51) 


Integration over a section through the flow field of a free 
jet in a plane perpendicular to the axis of flow leads to 


ae ee 
(pu) y Ag = Jo pu 2ur dr Gas) 
This merely states the conservation of the flux of ¥ and is 
an obvious consequence of neglecting the generation terms 


Sra ao peat 
By analogy to Eq. 49, it is assumed that 


We =- Aap (53) 


33 
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where 


Wee aks (54) 
4 OB toler 


and Oy is a function of x only. Equation 51 becomes 
Opu £ Ay oO on ) 
ae Oy. Aeipehe tele 


A solution of this equation is 


p (85) 
°y 
It has been observed'’°’ that successive profiles of ‘the 
momentum flux density ratio are similar. Assuming this 
relation to be general for the transport properties 
(pu) 
nase = fo 
(uly (56) 
From; iq...55 
(yu), r A Oy ea 
(bu), ‘ ras Y (57) 
But Eqs. 56 and 57 are compatible only if 
by = Cyx (58) 


Substituting Eq. 55 into 52 and integrating leads to 
Ky = (yu) 9 dg/ 


Therefore for a free jet issuing from a point source, 
Eq. 55 takes the form 


(Yu)o gg en (r/C 


; nC? x? 


x)? (59) 
y 


That Ls) Eq. 59 gives the distribution of the flux of any 
entity which originates wholly in the nozzle discharge 
based on the following assumptions: (1) the jet issues 


from a point source; (2) successive distribution profiles 
are Similar; (3) Eq. 53 is valid; 


10. Principle of Superposition 


Upon substituting Eqs. 53, 54, and 58 into Eq. 51 there 
is obtained 
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du Cy2x 0 dyu 


rah or 3p 


me (60) 


fiteeequation isslinear in Wu, and therefore a linear 
combination of particular solutions is also a solution. 
Thus, the distribution of yu in a flow field formed by the 
confluence of any number of point source jets issuing in 
parallel can be described by superimposing the solutions 
for the flow fields of the individual jets. 

For point sources discharging parallel to one another 


pu = : ee, (61) 


In this summation it is not necessary that all the point 
sources lie in the same plane perpendicular to the axis of 
flow provided that no point source impedes the flow from a 
source lying behind it. Furthermore, if a finite source may 
be taken as an infinite number of elemental point sources, 
then the summation may be replaced by the integration 
2, 
e~'rga/ Cutan) GA 


(pu) oy 
2 


2 
myx a4 (62) 


—— A 


“47. 0 


0 


Il. Distribution of Momentum Flux 


For the case of momentum flux, is pu. Substituting 
Bliis in Eq. 59 gives 


ee As Fey lrhc,x)? 


ato i Oa (63) 


The momentum flux density is measured approximately by a 
total-head impact tube, and Eq. 63 can be compared with 
such measurements. 

In cases of mass and heat transfer.at is frequently 
necessary to employ the mean velocity uw rather than the 
momentum flux. Introducing the mean and the fluctuating 
components into Eq. 63 yields 


Set 9 
(ot+o' )(utu' ) Ay 
Gout le 0 Oa ee 


? 


en (r/C,x)? 


Expanding the left-hand member, neglecting third-order 
moments, and recognizing that all terms containing only one 
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fluctuating component vanish, there results 


(po wt ys A 90" u’ ul? Ay ee )2 
—_—$—$$———— — — = = \ ry" x 
paws ce iene Se aCe é bead 


Following Corrsin,‘*°’ the density fluctuations are ex- 
pressed in terms of the temperature and concentration 
fluctuations. For an ideal gas in turbulent flow p = MP/RT 
to a first degree of approximation (see Section 7). 
Expanding Ap in Taylor’s series and neglecting terms 
of higher order 

MP P 


Ap =- 
zs RT 


For a two-component mixture 


mM =>~—1 
wey Lee 
u, +m, 
whence 
AM oes MAw . 
es Mo/M., 
i 1-M,/M, 
Substituting for AM in Eq. 65 and setting Ap, AT, and 
pw equal top", t, w' . respectively 
ed a f 
p stip ieee oes (66) 
Tw’ 5+G 
where 
MM, 
1-M)/M , 


Substitution of Fq. 66 into Eq. 64 gives rise to the terms 
t'u'/T u, w',u'/(w;, - G)u, and u!2/u2, The terms v' 2/y?2 


ye) 
and w! /u* also appear in slightly different circumstances, 
Corrsin discusses methods of measuring these by means 
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of a hot-wire anemometer. The mean temperature, T, can be 
measured by a thermocouple. The mean concentration of jet 
fluid, w,, can be obtained by standard methods of gas 
analysis. The mean density p can be calculated to a first 
@esree (of approximation from the mean concentration and 
temperature by the use of a suitable equation of state. 
Thus, methods of evaluating all quantities in Eq. 64 are 
available and u may be calculated from 


au 2AP ‘3 
ie [ |? 
ior hi ul? 
Olga os 
[evetin — 3 


where AP is the total-head impact tube measurement. 

Substituting md*/4 for Ag in Eq. 63 and setting r equal 
to zero gives the axial distribution of momentum flux for a 
circular nozzle as 


UG vs oe (67) 
ee 2C_,x 


Substitution of this into Eq. 63 gives 


2 
(ou ie a e-(r/c_ 2)" 


(ene) (68) 


which expresses the similarity of the flux profiles when 6 
Meobinear in x. 


12. Distribution of Flux of Sensible Heat 


For the flux density of sensible heat under conditions 
similar to those just discussed for momentum flux, ? is 
pC_t. The generation terms in Eq. 2 were neglected, and 
therefore changes in t due to chemical reaction are not 
considered. Likewise, phase transformations, Joule effects 
and the like are neglected. 

Substituting for Y in Eq. 59 leads to 


(oC tu) A 
ieee 8 8 Vc is 1? (69) 


5 2.2 
(pC, tu)y 9 mC) x 
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Here again, pC tu cannot be measured directly but, as 
before, one may write 


‘ re ‘ ‘ r , 9 
t 
(0Cpt nee PAG ET p Ay of t/Cpe) 


[1+ + + ] = 70 
(pc, tu) tcul 01ai eM 1 OP (70) 


where C, 1s assumed constant. Equation 66 applies to this 
case also. 


13. Distribution of Mass Flux 


In like manner for the mass flux density, ¥ is pw; 
and Eq. 59 becomes 


ORT ie Ay Poly foe 
(pw .u) 9 0 Go oar (71) 


Since there is no direct method of measuring pw.u, it is 
again necessary to introduce the mean and luctuseine 
components, giving the relation 


low, u), w" eu! piu’ w' Pp’ Ay it 7Gixed? 
———— [1+ — +-——t-——] = ———-e (72) 
(ow uy 9 wu pu w, p mG Sa 


Equation 66 again applies. 


1V. EXPERIMENTAL INVESTIGATION 
14. Apparatus 


The apparatus comprised a source of compressed air, a 
flow control and metering system, heat exchangers, flow 
nozzles for generating free jets, impact tubes, and a 
positioning mechanism for locating the impact tubes in the 
jets. A schemetic flow diagram is shown in Fig. 2. 

The compressed air was supplied by a diesel-powered two- 
stage air compressor having a capacity of 500 standard cu 
ft per min at a pressure of 100 lb per sq in. gage (Fig. 3). 
The compressor was connected to an aftercooler through a 


Moore Nullmatic 
Control Valve ——+O 


Vent to Measure 3 het or Cold Water 
Humidity 


Fisher Wizard 
Control! Valve 


ia From Compresser Orifice Plates 7 Flow Nozzles we 


><] 


Fig. 2. Schematic Diagram of Experimental Apparatus 


2-ft section of flexible hose which prevented vibrations 
from the compressor from traveling down the connecting 
lines to the experimental apparatus. The air passed 
through an oil and moisture separator into two air recelv- 
ers with a capacity of about 200 cu ft, and thence through 
two air filters consisting of porous ceramic discs 10 in. 
in diameter and 2? in. thick. 

The metering system is shown in Fig. 2. It consisted of 
nine interchangeable orifice plates covering a range of 
flow from 0.5 to 350 standard cu ft per min. The plate in 
use was inserted between two flanges fitted with three pins 
located asymmetrically so that the plate could be oriented 
only in one position. A bulb-type dial thermometer was 
used to measure the upstream temperature #2 deg F. The 
pressure on the upstream side of the orifice plate was 
measured with a Bourdon test gage giving readings accurate 
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Fig. 3. Air Compressor 


to % percent of the full scale, while the pressure on the 
downstream side was measured with a regular Bourdon gage. 
The orifice was operated in such a manner that the ratio of 
the downstream pressure to the upstream pressure was less 
than the critical pressure ratio for the flow of air under 
the conditions used. Hence the mass rate of discharge of 
air through the orifice was independent of small changes in 
downstream pressure, and the rate was determined entirely 
by the upstream pressure and temperature. 

For flow rates up to 50 standard cu ft per min the up- 
stream pressure was regulated by a Moore Nullmatic control 
valve; for rates of 50-350 standard cu ft per min, by a 
Fisher Wizard control valve (Fig. 4). Both valves covered 
the range of pressures from (0 to 100 lb per sq ln. gage. 

A part of the metered air from the orifice was passed 
through a battery of heat exchangers. This air was mixed 
with the main stream to adjust the temperature so that when 
it expanded adiabatically through the flow nozzle, the jet 
would be discharged at the ambient temperature of the room. 
Four shell-and-tube exchangers were connected in series, 
each having a transfer area of about 35 sq ft. These could 
be supplied either with hot boiler water at about 215 deg F 
or with cold water from the mains. Control of the main- 
stream temperature was accomplished by varying the fraction 
of the air sent through the exchangers. 
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Three manometers were used for measuring the impact 
pressures and the pressure immediately before the flow 
mozzle--a 40-in. mercury manometer, a 40-in. water. mano- 
meter, and a Wahlen gage.'’*) This latter gage wasused for 
measuring pressures from 0.001 to 1.000 in. of alcohol. 

The pressure control system, the heat exchanger, the 
flow control apparatus, and the instrument panel are 
shown in Fig. 4. 

Several flow nozzles having different throat diameters, 
throat lengths, and exterior designs were used in the 
experiments. All lad an ASME long-radius, low-ratio 
approach section. ‘*’ Two typical 0.898-in. nozzles with 
semi-streamlined exteriors are shown in Figs. 5 and 6. 

A mechanism for traversing in three dimensions was used 
to position the impact tubes in the flow field of the jets. 
This mechanism traversed the region within 100 cm axially 
from the nozzle, +30 cm taterally, and +25 cm vertically. 
It consisted of two rails having longitudinal V-notches in 
which a cross rail rested on ball bearing end-supports. 
The cross rail was surface-ground on two adjacent sides and 
provided a bearing surface for a cross carriage. Alignment 
of the V-notched rails with the nozzle axis was held within 


an angle of 0.2 deg. 


Fig. 4 Flow Control and Metering System 
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Fig. 5. Cross-section of the 0.898-in. Semi-Streamlined Nozzle 


The impact tubes, discussed in detail in Section 16, 
were mounted on a vertical rod, 7/8 in. in diameter, which 
passed through a sleeve clamp. The elevation of the impact 
tube was determined by a vertical cathetometer, while the 
axial and lateral coordinates were indicated by pointers 
bearing on centimeter scales attached to the V-notched 
rails and the cross rail, respectively. 


15. Calibration of Instruments 


The standard for the velocity and momentum measurements 
consisted of two orifice plates which were standardized in 
the laboratories of the manufacturer. Water, measured 
volumetrically to 1 part in 10,000 in carefully calibrated 
tanks, was passed through the orifices at constant rates 
for periods of time measured to 1] part in 3000. The 
pressure difference across the orifices was determined as a 
function of the flow rate in the. same range of Reynolds 
Number later used with air, 

The orificesiwere ‘then used’ to. cali bratesaetarnieme 
steel ASME standard flow nozzle having a throat diameter of 
0.952 in. The nozzle was calibrated by connecting it in 
series with the standardized orifice and an air source. 
The, pressureydifierence across. the orafpemendvacrocs the 
nozzle, along with the corresponding temperatures, were 
observed at various flow rates; the discharge capacity of 
the nozzle was expressed in terms of the pressure and 


temperature just before the converging section by means of 
the following equation '‘?’ 
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h P 
Oeers * se Maal tT (73) 


By a series of 42 calibration experiments in which the 
pressure differences across the nozzle were less than the 
critical difference, the following relation was found 


Q = 3 G@osserelye °}* rd 7 
aio oe See ae ce (74) 


which covered air flow rates from 0.06 to 5 standard cu ft 
per sec. In another series of 12 experiments in which the 
pressure differences were above the critical difference, 
the relation 


P 
4 POA Cereal | 75 
Qscfs = 3.0286 (Sagar (75) 


was found covering the flow rates from 5 to 8 standard cu 
ft per sec. The standard flow nozzle was used to calibrate 
the orifice meter described in the previous section, 
Calibration of the impact tubes is discussed in the next 
section. 


q 


Figs 6. Exterior View of the 0.898-in. Nozzle 
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16. Total-Head Impact Tube 


a. Calibration 

The total-head impact tubes were made from No. 16 
hypodermic-needle tubing having an outside diameter of 
0.065 in. and an inside diameter of 0.050 in. Before use 
the ends were carefully squared and burrs were removed. 
In addition the tubes were calibrated and their readings 
interpreted according to the equation 


2 
AP = = oe 
epee (74) 


From reference (3) C, may differ from unity and should be 
determined experimentally. Goldstein‘?°’ defines a cali- 
bration coefficient for pitot tubes by the equation, 
AP = Kou?/2. Hence, K = 1/C,2. He reports that the value 
of K at speeds varying from 20 to 60 ft per sec did not 
vary from unity by more than +0.1 percent, while the value 
at speeds from 6 to 20 ft per sec did not differ from unity 
by more than +] percent. 

There are two explanations for values of C, below unity. 
Barker'*’ points out that at low velocities impact tubes 
of small radius are influenced by a viscosity effect, which 
may be taken into account by the equation 


where u is the velocity, yu the gas viscosity, and r the 
internal radius of the impact tube. It follows that the 
coefficient Cy at low velocities is given by 


pur 
2 
2 
C2 =— 
POs eo a ae 
2 Dior 


Using this equation with r = .050 and w for air at 70 deg F 
and atmospheric pressure, the values for C, were calcula- 
ted as a function of u and are given below. Impact tube 
coefficients of this order of magnitude were also reported 
by Homann'3°’ under similar conditions. 
At high velocities the stagnation pressure, static 
pressure, and kinetic head are related by the equation!*?) 
P,-P ] 1 
= 1 +—M? +—y4', : 


vag 
q A 40 Teoden 


a 
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Values of Calibration Coefficient 
at Low Velocities 


u 
ft/sec Cy 
4 0.971 
10 0.988 
20 0.994. 
50 0.998 


where P, is the stagnation pressure, P the static pressure, 


q the dynamic pressure, and M the Mach number. The quan- 
tity (P,) - P) is the pressure difference AP indicated by an 
impact tube, and the dynamic pressure q is given by pu‘“/2. 
Therefore 

20P ae Sere! = fab ne. 1 MOgeee 
pu oF q 4 40 #1600 


Using this relation C, can be calculated as a function of 
u, and the values so Shrained are given in the table below. 


Values of Calibration Coefficient 
for High Velocities 


u 

ft/sec GF 
100 0.999 
200 0.995 
300 0.990 
400 0.983 
500 0.975 
600 0.963 
700 0.949 
800 0.934 
900 0.918 
1000 0.901 


These values of C, for both low and high velocities are 
shown in Fig. 7 as a function of the velocities. 


6. Effect of Angle of Attack 


The "ideal" impact tube is defined as an instrument 
which measures the pressure at a point of stagnation in a 
flowing stream. Since the velocity is zero at such a point, 
the ideal impact tube has characteristics whichare indepen- 
dent of the direction it faces; that is, the pressure at a 
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point of stagnation is independent of the orientation of 
the obstacle producing the stagnation. 
Actual impact tubes, however, do have directional 
characteristics. This fact was demonstrated in an experi- 
ment in which the axis of the tube was rotated through 
an angle of 90 deg with respect to the direction of flow 
of a nonturbulent stream. Up to an angle of about 15 deg 


1.0 


OSES 


Cr 


0.90 


1 5 10 50 100 500 1000 
Velocity, Feet per Second 


Fig. 7. Total-Head Impact Tube Coefficients 


the indicated pressure was substantially constant. Beyond 
that angle the pressure decreased, becoming less than 
atmospheric for angles greater than 60 deg. 

Since in a turbulent stream the direction of flow at 
the impact tube fluctuates, the question arises as to 
whether the average angle of attack is large enough to have 
an effect on the pressure readings. The ratio of the 
root-mean-square of the transverse velocity fluctuations to 


root-mean-square of the axial velocity, V v'2/y2, may be 
regarded as roughly proportional to the sine of the 
root-mean-square of the average angle of attack on the 
impact tube. Corrsin'*°’ found this ratio to be nearly 0,2 
across the entire flow field of a free jet, corresponding 
to an average angle of attack of about 12 deg. 

Since this average angle is less than 15 deg, it was 
concluded that impact tubes used in free jets are not 
subject to significant errors which are due to fluctuations 
in the angle of attack. 
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ce. Effect of Turbulence 


In actual use the momentum flux densities, pV2, were 
measured throughout the flow field of turbulent jets by 
“impact tubes, and the calibration coefficient defined by 
Eq. 76 was applied regardless of the level of turbulence 
and hence regardless of the velocity fluctuations ex1sting 
at the point of measurement. Other investigators have 
followed essentially the same procedure, because the basis 
for taking the effect of the velocity fluctuations into 
account is not well understood. In an attempt tc determine 
the effect of the velocity fluctuations, the following 
analysis and experiments were performed. 

For steady-state flow, changes in static pressure in any 
direction may be represented in vector notation by 


dP = YP+dR Cra) 


where R is a vector distance. When the momentum flux 
Variables are substituted into Eq. 2 and the first two 
terms for unsteady flow and viscous shear are neglected, 
there results 


-V+(pV:S)V-VP*S = 0 


Successive expansion in the x, y, and z directions and 
recombination lead to an expression for VP. Substituting 
‘this expression for VP into Eq. 77 yields 


dP =(V:ouVi+V*pvVj+V*owVk) *dR 
Finally application of the equation of continuity gives 

dP = (oVVsuitoVV:vj+oVV-wk) dR 

Consider now an instantaneous streamline through a 

stagnation point and assume that the distance along the 
streamline from a point where the velocity has the free- 
stream value to the stagnation point is so short that this 
streamline is substantially a straight line between these 
points. If the x-axis is oriented parallel to this stag- 
nation streamline, then the velocity components normal to 
it, v, and »,, are zero and 


d 
dP = pu, “& dr (80) 


ax 
8 


where both the distance, x,, and the velocity, u,, are 
measured along the stagnation streamline. Sinecesu 38 
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equal to the magnitude of the velocity vector, integration 
of Eq. 80 along the streamline from the free-stream point 
to the stagnation point, when averaged with respect to 
time, gives 

AP = doy? 

2 (81) 

This expression may be expanded in terms of the mean and 
fluctuating components to give, for incompressible flow, 


BP = S-plu2tu!24y!24u'?) (82) 


and solving for u yields 


QAP eller to et ee ae 
a = ul tay! 2-2) (83) 


Unset 


Hence strictly speaking, u may be obtained from impact tube 
measurements in turbulent streams only if the magnitudes 
of the fluctuating components are known. 

In addition, the momentum flux density in the x direction 
1s given by 


pu2 = 2AP ~ pv'? - pw!? (@4) 


and substitution of 0 from Eq. 81 into this equation 
results in 


oe —Agte er omens 
pus = QAP(l-—--——)} (85) 
v2 y2 


Data concerning w'2/y2 are_not available; however, if this 
1s approximately equal to u'2/V2, then the term in paren- 
theses may be estimated from Corrsin’s‘?°! work to be 0.92 
along the axis of a free jet. From this analysis it would 
appear, therefore, that in determining the x-directed 
momentum flux, the correction to be applied to the pressure- 
difference reading on the axis, in order to allow for the 
velocity fluctuations in the y and z directions caused by 
turbulence, amounts to about 8 percent of the reading. 

As a more direct means of determining the effect of 
turbulence on the impact tube measurements, several other 
experiments were performed in which impact tube measure- 
ments in nonturbulent and varying levels of turbulent 
jets were compared, 

The nozzle used in these experiments is shown in Fig. 8. 
It had a throat diameter of 0.75 in., and the converging 
section followed ASME specifications for long-radius low 
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Fig. 8. Cross-section of the 0.75-in. Nozzle 
with Screen and 0.25-in. Throat Extension 


ratio flow nozzles.‘%’ The nozzle was equipped with four 
removable throat-extension sections -- 0.25, 0.50, 0.75, 
and 1.00 in. long. For generating a turbulent jet a 


35-mesh copper screen having a wire diameter of 0.243 mm, 
Square openings of 0,44] + 0.030 mm on a side, and a free- 
flow area of 41.3 percent was inserted at the entrance to 
the throat section. To determine the scale and the level 
of turbulence at the discharge end of the throat section, 
reference was made to empirical equations developed from 
similar work by Dryden. ‘**’ These values are tabulated 
below as a function of the length of the throat section. 

In the experiments themselves the impact tube was loca- 
ted at the discharge of each throat section in turn while 
the air flow rate through the nozzle was maintained at the 
constant value of 180 ft per sec. To obtain measurements 
in a nonturbulent jet for comparison, the experiment with 
the 0.25-in. throat section was repeated with the screen 
removed. Although the results are not presented here, they 
were in fairly good agreement with the qualitative predic- 
tions from Eq. 82. With the mean velocity held constant 


Turbulence in Discharge from a Nozzle with Screen 


Length of section Turbulence level Scale of turbulence 
in. percent mm 
0.25 uo 0.12 
0.50 4.2 O13 
0.79 ou Oieko 
1.00 225 0.16 
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while the turbulence level is raised, there is a definite 
effect on the impact tube readings, but this may be partly 
caused by minor changes in the velocity distribution 
across the nozzle. 

Consequently in accordance with the common practice of 
other investigators, the impact tube measurements were not 
corrected for the fluctuating components of velocity 
because the basis for making proper corrections is not well 
understood. Although these fluctuations in velocity may be 
one of the chief sources of error, it is believed that this 
error 1S not exceSsive in most regions of the jet. To 
correct the impact tube measurements, therefore, the 
calibration coefficient from Fq. 76 was used. 


17. Single Free Jets 


The experiments described in this section consisted of 
obtaining impact tube traverses along the axis of single 
jets and also along a radius of several cross sections of 
the jet at various distances from the nozzle. Several 
nozzles having different diameters were used including the 
special nozzle described in the previous section. The 
velocity at the nozzle was also varied in the experiments. 


a. Transport of Momentum Flux 


The results of several measurements of the momentum flux 
density in a free jet of air from a 0.898-in. nozzle are 
shown in Fig. 9. Two sets of data are presented in the 
figure. One set (Table 1) was obtained by making radial 
impact tube traverses at axial distances of 10, 15); 320s 
and 30 nozzle diameters while the discharge velocity was 
held constant at 403.1 ft per sec. The other set (Table 2) 
was obtained by making a radial traverse at 20 nozzle 
diameters at each of four discharge velocities ranging from 


177.9 to 803.6 ft per sec. 


In Fig. 9 these data are compared with the probability 
distribution of momentum flux which was used to arrive at 
Reichardt’s hypothesis. For this comparison Eq. 57 was 
modified as follows. By defining ry as the radial distance 
at_which (fu), is ‘one-half’ the value oh thesaxae 
(tu), fuetoliows that 


ef ry/ by)? Bi. 
ys 
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Fig. 9. Comparison of Radial Profiles of Momentum Flux Density 
Ratio in a Single Free Jet with the Probability Distribution 


from which 


2 
ry 
b 2 = 72 
z Lae2 
Substituting for by? in Eq. 57 gives 
(Gale, s eo (ry)? ine (86) 
(Yu), 9 


Thus Eq. 86 implies that the ration(¥u), /(¥u), 4 1s 


unique function of r/r, and that this function contains no 


arbitrary constants regardless of the functional dependence 


of by on x. 


For the case of momentum flux ¥ is pu, therefore Eq. 86 


becomes 


Bere, M eo(ty)21n2 (87) 


a3 
(pu V6 


5| 
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Calculated values of Eq. 87 are represented by the solid 
line in Fig. 9. This line, the probability curve, corre- 
lates the data reasonably well. The points have a slight 
tendency to lie below the curve in the range of r/ry from 
0 to 1.0 and to lie above it elsewhere. 


6. Effect of Similarity Conditions on Distribution of Momentum Flux 


In Section 9 the similarity of radial profiles was 
assumed and was used to deduce the relation between 6 and 
x. There are three kinds of similarity involved in Fq. 63. 
It requires that successive profiles of momentum flux must 
be similar (geometric similarity), and that the ratio 


(pu?) PARTE ean must be independent of the discharge 


velocity (dynamic similarity) and also of the nozzle 
diameter (dimensional similarity). 

The results of several impact tube measurements taken 
through a free jet issuing from a nozzle at a velocity of 
403.1 ft per sec are shown in Table 3 and Fig. 10. As 
suggested by Eq. 59, the negative logarithm of the momentum 
flux density ratio is plotted against (r/x)? for traverses 
taken at various distances from the nozzle. 
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Fig. 10. Geometric Similarity: Variation of Negative Logarithm of 
Axial Momentum Flux Density Ratio with Square of Radial Distance Ratio 
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The linearity of the individual sets of data is reason- 
ably good. However there is some trend in the slopes with 
distance from the nozzle. Values of the spreading coeffi- 
cient, C,, calculated from these curves are tabulated 
below. The mean value of C, is about 0.075, and the 
maximum deviation is about 4 percent. The values at first 
decrease and then increase with increasing distance from 
the nozzle. The variation is not large, however, and the 
assumption of geometric similarity is good. 


Axial distance ratio Momentum spreading coefficient 
x/d Ge 
10 0.078 
15 0.074 
20 0.073 
25 0.074 
30 0.076 


* 


In another series of experiments the impact tube was 
placed at a fixed point on the axis of flow while the dis- 
charge velocity was varied. Values of the ratio of the 
momentum flux density on the axis to that at the nozzle‘)?! 
are presented in Table 4 and are plotted against discharge 
MeLOcr GY in, Fig: kl. Lhe change in the momentum flux 
density ratio with velocity is small only at velocities 
in excess of 200 ft per sec and at distances from the 


nozzle greater than 20 nozzle diameters. At efflux ve- 
locities above 200 ft per sec the variation in C1, as 
calculated from these curves, is small. The values lie 


between 0.070 and 0.075. 


Table 4 
Dynamic Similarity: Variation of Axial Momentum Flux Density Ratio 
with Velocity at Nozzle 
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Density Ratio with Velocity of Discharge at Nozzle 


Finally the effect of varying the nozzle diameter on the 
momentum flux density ratio at a constant value of x/d is 
given in Table 5 and shown in Fig...12. Although the scat- 
ter is large, it would appear that there is no particular 
trend in the ratio with nozzle diameter. 

It was concluded from the foregoing series of experi- 
ments that reasonably good geometric and dimensional 
similarity exists in the free jets at a given velocity, but 
that comparisons between jets of different discharge velo- 
cities should be made with caution. Actually the existence 
of similarity is not necessary because the application of 
Reichardt’s hypothesis does not in any way depend on the 
conditions of similarity. If geometric and dimensional 
similarity do not exist however, the analysis of Section 9 
would be considerably more complicated. 


c. Effect of Induced Turbulence on Spread of Jets 


To determine whether induced turbulence has any effect 
on the dispersion of a jet, a series of experiments were 
performed using the special nozzle (Section 16) with the 
turbulence generating screen in its throat section, The 
turbulence generated by the screen was presumed to persist 
in the normally nonturbulent core, and was calculated by 
Dryden’s method‘**’ to be approximately 8 percent. The 
0.25-in. throat extension was used for al] experiments. 
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6 and 20. 
nozzle. 


versus r/d. 
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Table 5 
Dimensional Similarity: Variation of Axial Momentum Flux 
Density Ratio with Diameter of Nozzle 


ft/sec 


226.76 
219.308 
218.50 
198.35 


0.051 


The data were obtained by making traverses with the 
total-head impact tube at axial distance ratios, x/d, Of 2s 
For comparison purposes measurements were also 
made with the turbulence generating screen removed from the 
Both these sets cf data are shown in Table 6 and 
are plotted in Fig. 13 as the momentum flux density ratio 
The nonturbulent measurements are shown as the 


solid curve. 
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Fig. 12. Dimensional Similarity: Variation of Axial Momentum 


Flux Density Ratio with Nozzle Diameter 
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From Fig. 13 it is apparent that the agreement at 
2 nozzle-diameters is very poor while at 6 nozzle diameters 
and Beyond ibvisiquite good: “This does not» mean that the 
induced turbulence is causing the jet to spread, however; 
“because if that were the case, then it would be expected 
that the discrepancy at 2 nozzle diameters would also be 
present at 6 and 20 nozzle diameters. It would appear 
therefore that the disagreement at 2 nozzle diameters is a 
result of the effect of the turbulence on the impact tube 
rather than an indication of the spread of the jet. 


d. Transport of Sensible Heat Flux 


A criterion for recognizing situations in which the 
transport of heat, mass and momentum is similar is provided 
by Eqs. 1 and 2. Similarity exists when the same terms ap- 
pear in these equations for the different kinds of trans- 
port. Two examples will illustrate different cases. 

First, consider the steady-state flow of an incompres- 
sible fluid in a pipe. The term for accumulation drops out 
and the two terms for transport by molecular processes and 
by fluid motion are recognized as similar for all three 
cases. There is a difference, however, in the generation 
term. Ordinarily there is no generation of heat or mass in 
the tube, and this term is dropped. However, the pressure 
drop along the pipe results in the generation of momentum 
equal to that transferred to the pipe wall, and this term 
1s required when momentum is considered. Because of the 
inclusion of this generation term in the momentum equation 
and its absence from the heat and mass transport equations, 
these three transport equations do not lead to similar 
solutions. A physical situation can be imagined, however, 
in which they do. Thus, if a fine metallic powder were 
suspended in the fluid and heat were generated by inductive 
heating uniformly throughout the mass of the fluid and 
at a rate equal to the rate of loss of heat to the tube 
walls, then the generation term required in the heat 
transfer equation would be similar to that for the transfer 
of momentum. 

As a second case consider the steady-state flow of an 
incompressible fluid in a free jet. The terms for accumu- 
lation and transport by molecular processes and fluid 
motion are again handled as for the flow in a tube. How- 
ever, the pressure throughout a free jet 1s uniform and 
hence there is no provision for generating momentum. The 
generation term is now omitted from the momentum equation 
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Table 6 
Radial Profiles of Axial Momentum Flux Density Ratio 
in a Free Jet Discharged Through Nozzle with and without Screen 
d = 0.75 in.; average velocity at nozzle = 180 ft/sec 
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Table 6 (Concluded) 
Radial Profiles of Axial Momentum Flux Density Ratio 
in a Free Jet Discharged Through Nozzle with and without Screen 
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as well as from the equations for heat and mass transport, 
and all three transport situations are similar. 

No measurements were made in this work with a tempera- 
ture gradient in the jet. However, Hinze‘’?’ has reported 
extensive measurements on the temperature and momentum flux 
distributions in a jet of air 30 deg C above and below the 
ambient temperature. From his impact tube measurements he 
calculated u by a method similar to that described in 
Section ll, using data obtained by Corrsin. ‘7°? Hinze 
made a large number of measurements, and in these there is 
relatively good agreement between the distribution of 
temperature and momentum flux with only moderate dispersion. 

While Eq. 63 can be compared directly with experimental 
data, this is not true of Eq. 69. The reason 1s of course 
that the pitot tube measures, as a first approximation, the 
total momentum flux; whereas there 1s no instrument which 
measures the total heat flux. A thermocouple measures only 
the mean temperature, not the fluctuations. A better appre- 
ciation of the problem of making measurements 1s obtained 
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Fige 13. Radial Profiles of Momentum Flux Density Ratio in a Free 
Jet Discharged through a Nozzle with and without Screen 


by inspecting Eq. 70, which is an expansion of Eq. 69. 
A complete comparison between experiment and theory could 
be obtained if the turbulent parameters in the bracket 
of Eq. 70 could be measured. Such measurements have, 
in fact, been made by Corrsin by means of the hot-wire 
anemometer. However, the points for t’u’ scatter so widely 
that the data cannot be used for the present purpose. 

An approximate comparison of Eq. 70 with the experimen- 
tal measurements can nevertheless be made by recognizing 
that the value of the quantities in the brackets of Eqs. 64 
and 70 lies between ] and 1,2. These equations may be com- 
bined by first dividing Eq. 70 at x,r_by the same equation 
at x,o and then eliminating oer and Ue 0 by using Eq. 64. 
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In these measurements the range of temperature was small 
and the density p essentially uniform. The resulting equa- 
tion contains the ratio of the bracketed quantity in ive rete 
to the square root of the bracketed quantity in Fq.(64. As 


a first approximation, this ratio may be assumed to be 
unity. The resulting equation 


t ih il ro? 
pokes (ee ees 

<5 2 

aa 2C,, C,? ft 


then gives the distribution of the mean temperature as a 
function of position and does not involve turbulent quanti- 
ties. Hence, this equation may now be compared with 
Hinze’s measurements. The temperature ratio appearing on 
the left-hand side of the equation when pee versus the 
square of the dimensionless radius (r/x)“ on semilog paper 
should be a straight line. Figure 14 shows such a plot. 
The agreement is clearly excellent, indicating that the 
transport of heat may be treated like that of momentum, 
The dispersion of the large number of measurements is shown 
by the shaded area, while the straight line is a good mean 
over most of the range. 

In a jet of hot fluid there is a radial flux of heat 
‘>roportional to t‘v’. If the temperature and velocity 
luctuations were independent of one another, this quantity 
vould vanish and the turbulence would not be conducive to 
-he radial heat flux. Actually, these fluctuations are not 
“independent. A velocity fluctuation at any point may be 
“sisualized as a consequence of the arrival of an eddy at 
‘that point. The temperature of this eddy is representative 

of the average temperature at its starting point. As a 
result the velocity fluctuation is accompanied by a tempera- 
‘ture fluctuation. Furthermore, if the velocity and tempera- 
ture gradients have the same sign (es  inta free jet with 
the jet at a temperature higher than the surroundings), the 
fluctuations of these quantities should also have the same 
sign. Similar considerations apply to the axial flux t’u’ 
appearing in Eq. 70. 

On the other hand, since the density 1s inversely 
proportional to the temperature, the most probable combina ~ 
tions of density and velocity (or temperature) fluctuations 
are those in which the fluctuations have opposite sign, and 
the terms p'u’ and p’t’ are negative. It is possible that 
the sum of these terms could be greater in magnitude than 
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Fige 14. Heat and Mass Flux Density Ratio Profiles 
in a Free Jet, from Hinze's Measurements 


t’u’ in which case the quantity enclosed in brackets in 


Eq. 70 would be less than unity. Comparing Eqs. 69 and 70, 
it is seen that 


(p: C. 4 Le 1] 
5 eee 
TR ae Pee 
Gat. 1, 


tu pu pt 


A value of less than unity for the quantity in brackets 
means@that the heat flux density calculated from the mean 
components of flow, pC,t u, is greater than the mean value 
of the instantaneous heat flux density, pC tu. In the case 
of sufficiently, cold jets these consmteretvons would be 
reversed, and the flux density calculated from the mean 
components would be less than the true mean value. 
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_ For the case of a heated jet where the temperature 
difference with respect to the ambient condition is small 
the terms p'u' and p't' are probably smallcompared to the 
term t’u’ inasmuch as p is inversely proportional to the 
absolute temperature, Thus the quantity in brackets should 
exceed unity, and pC,t u should be less than pC tu. 
The order of magnitude of the fluctuating terms can be 
estimated by integrating measured values of pC_t u over 
the flow field in a plane perpendicular to the axis of 
flow and comparing this total flux of heat to that flowing 
out of the nozzle. 


The integral 


dr (88) 


(pC, tu)y 9 Ao 


was evaluated from Hinze’s data, and the most probable 
value was found to be 0.73; 27 percent of the total heat 
flux, then, was not accounted for in the measurement of 
the mean quantities and was transported by the fluctuating 
components. Hence the determination of heat flux re- 
quires that both the mean and fluctuating components 
be taken into account. 


e. Transport of Mass Flux 


The criterion for the similarity between the transport 
of mass and momentum and between mass and heat is like that 
for the similarity between the transport of heat and momen- 
tum described in Section 17d. This criterion applied to a 
free jet makes it evident that the transport of mass should 
follow a pattern like that for the transport of momentum. 
Such a case is that of a nozzle discharging a fluid which 
is less dense than the surrounding fluid. Under these 
conditions Eq. 72 is applicable. 

Although no experiments were performed with nonhomogen- 
eous jets, Hinze '??) has reported extensive measurements 
of concentration and momentum flux density profiles in a 
jet of air containing 1.1 percent of "town gas." These 
data may be compared with Eq. 72 by a procedure which par- 
allels that used to compare the temperature distribution 

Ww 


with Eq. 70. In this case In—tt would be expected to be 
Ws. 
0 
linear in (aaj). This concentration ratio was found by 


Hinze to fall very nearly in the same range as the tempera- 


ture ratio. In fact the same curve shown in Fig. 14 may 
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be used for both ratios withequally excellent agreement. 
The comparison indicates that the transport of mass may be 
treated like that of momentum or heat. 


An equation entirely parallel to Eq. 88 may be written 
for this case 


(po w, u) 2rr 


~0 u tat 


The profiles of the product of the mean concentration and 
mean velocity were fitted with an error curve, and the in- 
tegration of this equation was performed analytically. The 
value of F, was estimated to be 0.79 indicating that 2] 
percent of the total mass flux can be attributed to the 
fluctuating components. 

The jet fluid used by Hinze had a density negligibly 
different from that of air so that the terms in Eq.%@ 
involving o' were negligible. The quantity in brackets 
should have exceeded unity, as was found experimentally. 
Values less than unity might be expected from jets of 
hydrogen or helium. 

The experimental evidence described above indicates that 
the Reichardt hypothesis as extended in Fq. 53 to the 
general case of the transport of heat and mass as well as 
momentum applies to a simple free jet within the accuracy 
of the measurements. The practical utility of this prin- 
ciple would be greatly increased if it were shown to apply 
to more complex situations. Of interest are multiple jets, 
jets adjacent to solid surfaces, and jets issuing into 
tubes. A still more complicated case is that of a free 
jet in which there is a simultaneous transport of heat, 
mass and momentum, accompanied by a chemical reaction 
(as in a flame) or by the evaporation of droplets (as 
in spray drying). Several such situations are discussed 
briefly in Section 19. 


18. Application of the Principle of Superposition to Jets 
a. Single Free Jets 


Total-head impact tube measurements were taken in the 
flow field of a free jet in the region adjacent to a nozzle 
which was designed according to ASME specifications and 
with an inside diameter of 0.898 in. Air was discharged 
at a velocity of 362 ft per sec at the temperature of 
513 deg P. The results of these measurements are shown 
in Fig. 15 and Table 7. 
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Figs 15. Radial Profiles of Momentum Flux Density Ratio 
Adjacent to a Finite Nozzle 
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Equation 59 applies to a jet issuing from a point 
source. The correlation provided by Eq. 63 fails in the 
region adjacent to the nozzle--approximately closer than 
10 nozzle diameters. However as a consequence of Jeich- 
ardt’s hypothesis (Eq. 49), the equation of motion in the 
axial direction (Eq. 48) is linearized, and therefore a sum 


Table 7 
Radial Profiles of Momentum Flux Density Ratio 
Adjacent to a Finite Nozzle 
d = 0.898 in.; uy = 363 ft/sec 
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RBPrRreOoOcCcococoacacoccocoocooaeo 
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$96 
000 0.268 0.982 0.981 
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000 0.373 0.835 0.836 
0.403 0.756 0.756 
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0.610 0.044 0.0451 
0. 623 0.027 0.0275 
0.636 0.015 0.0153 
0.645 0.009 0.00913 
0.654 0. 0027 0.00273 
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of particular solutions (Eq. 48) 1s also a solution. 
This principle of superposition may be applied to com- 
binations of point source jets provided that the boundary 
conditions can be met and that there are no nonlinear 
interactions among the jets. 

For the case of momentum flux, 2 1s pu and Eq. 62 may be 
written as 


(pu?) Spans 
(ou?) sro dA | r gal CnX ‘i 
=! 0 mC 2,2 (89) 


In applying this equation to the flow adjacent to a 
finite nozzle a question arises concerning the nonturbulent 
cone described in Section 3 and shown in Fig. 1. Peich- 
ardt’s hypothesis is obviously valid where there is no 
velocity gradient perpendicular to the mean flow, as is the 
case inside this nonturbulent cone. It may not hold 
however for the transition region at the boundary of this 
cone. However, by neglecting this possibility, the plane 
of discharge of a finite nozzle may be treated as comprising 
an infinite number of elemental point source jets. The sum 
of the momentum flux from all of these at a point (x,r) is 


A ( 2 eotets zee cosy 
Su) a LaF pu OF. 3 CG x 
EDs 70 EP BA rere aac : a ee 
7™Cm Xx 


(90) 


where (s,y) is a system of polar coordinates in the plane 
of the nozzle discharge. For the case where (pu 0°. tae 
constant, integration gives : 


r d r 2n d 2p 
2 “Uaioe ISS laee ee) fee 
(ou hoe ieee Cx + IC x! ] “ 3G) : Cx) 
(pu )y 4 n=O n! es p! 


The series converges slowly; when r = x = d/2, the 45th 
term is the largest, and over 100 terms would have to be 


evaluated. However, a simple expression for the axial 
velocity may be obtained by setting r equal to zero 
(pu2) ={ : ? 
fie Ole, Lo acewnl eeiae (91) 
(ou?) 4 4 
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At distances from th 
e nozzle such that d/2C << i 
reduces to foals era 


(pu*), 0 d 9 
pan aS IT 
ou 1s, 0 


Which is the well-known parabolic law for the distribution 
of momentum along the axis of a free jet in the zone of 
fully developed turbulence. '2+19+58, 8,70) 

: It 1s more convenient to evaluate Eq. 90 by graphical 
methods than to use the integrated equation. For this 
‘purpose the origin of coordinates (s,y) is shifted to the 
point (x,r). Then when (pu2), 2 is constant over s, 


Eq. 90 becomes 


2 sharp (5 a)? 
(pu*) . ae dpe \2 elc,x/d) 
tii mae 
(pu?) , 5 d 2 
_y,-(r/d)?2 + (s/d)?-% 2 
EGE eee 
cos“ f9(;/d)(s/d) aay (92) 


‘A single graphical integration suffices to evaluate the 
‘momentum flux density ratio. The integration was carried 
out for various values of r/d at values Ofer /ar OLel gos 4 OF 
8, and 10, with an assigned value of the spreading coeffi- 
Peient, C,, of 0.075. The results of‘the calculations are 
‘given in Table 8 and are shown as solid lines in Fig. 15. 
Tt is seen that reasonably good agreement between pre- 
diction and experiment was obtained, At axial distance 
ratios of four and greater, the momentum flux density 
ratio on the axis is greater than predicted. Elsewhere 
it tends to be lower. 

An estimate of the precision of the data may be obtained 
by integrating the experimentally determined distribution 
of momentum flux over the flow field, and dividing by the 
rate of flow of momentum out of the nozzle. Clearly, 
regardless of the relation between pu’ andr, the value of 
this ratio must be unity; ‘thus 


(our Q7r 


dre= 1 
(pu7) 99 Ao Z (93) 
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Table 8 
Values of Momentum Flux Density Ratio 
Calculated by Graphical Integration of Eq. 92 
(C,, =' 0,075) 


x/d 


RPrRocCcoococeo 
YANUOUANMWNSO 
ooooocdo 


Because of the assumptions involved in deriving Eq. 92 it 
satisfies Eq. 93. The experimentally determined distribu- 


tion of (Coan?) , must also satisfy Eq. 93 within the 
limits of experimental error, and the deviation of the 
ratio from unity is a measure of the reliability of the 
measurements. The integration indicated in Eq. 93 was 
performed by the use of Simpson’s rule; the results are in 
the table below. The maximum deviation is of the order of 


4 percent, and there is a slight trend from higher to lower 
values with increasing x/d. 


Axial Distance Experimental Momentum Flux Ratio, 
Ratio, x/d Integral 


1 1.04 
1.04 
1.03 
1.02 
0.99 
0.98 


COON & NY 


— 


Bul. 4/3. MOMENTUM, MASSAAND HEAT IN TURBULENT JETS 73 


The discrepancy between experimental and predicted 
values of the momentum flux density ratio is apparently 
greatest along the axis of flow as shown in Figo 16. “Thus 
the agreement with Eq. 91 is poorer than with the pene 
correlation of the radial distribution, Eq. 87. Equation 
91 1s shown in Fig. 16 as a solid line. The experimental 
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Fige 16. Axial Profile of Momentum Flux Density Ratio 
Adjacent to a Finite Nozzle 


points fall on the line at values of x/d greater than 15 
and less than 3.5, but lie above the line elsewhere. The 
experimental values are given in Table 9. 


6. Two Parallel Jets 


The momentum flux density was measured in the flow field 
of two parallel jets. The nozzles are shown in Fig. 17. 
These had inside diameters of 0.898 in. and were located on 
3-in. centers in a plane perpendicular to the axis of flow. 
Air was discharged at a velocity of 388 ft per sec at am- 
bient temperature. Several traverses on a radius parallel 
to the line of-centers of the nozzles (¢ = 0°) are ‘shown 
in Fig. 18 and in Table 10. 

The principle of superposition may be applied to this 
-case of two parallel jets issuing at the same velocity from 
nozzles having the same internal diameter and lying in the 
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Table 9 
Axial Profile of Momentum Flux Density Ratio Adjacent to a Finite Nozzle 
d = 0.898 in., uy = 363 ft/sec, r/d = 0 


A p (pu), 9 
x/d AP. Das: 
0 2 

(pu 0,0 

0.67 1.000 1.000 
1228 1.000 1.000 
Tear 1.000 1.000 
PA sik 1.000 1.000 
2.54 0.999 0.999 
2.98 0.999 0.999 
3.42 0.994 0.994 
3.86 0.987 0.987 
4,30 0.972 0.972 
4.74 0.941 0.940 
5.17 0.904 0.906 
.61 0.856 0.858 
6.05 0.800 0.803 
6.49 0.745 0.750 
6.93 0.693 0.698 
7.36 0.641 0.646 
7.80 0.590 0.597 
8.25 0.542 0.549 
8.68 0.499 0.506 
9.56 0.431 0.438 
10.45 0.368 0,375 
11.34 0.318 0.325 
12.19 0.278 0.284 
13.08 0.248 0.253 
15.32 0.184 0.188 
17.54 0.145 0.149 
19.68 0.109 Oo112 
Pa eas 0. 087 0.089 


same plane perpendicular to the axis of flow. Equation 6] 
1s readily developed in simple functions. The origin of 
cylindrical coordinates (x,r,¢) is placed midway between 
the two nozzles, and the distance between nozzles is 
designated: ds D... Then for the momentum flux density 
ratio Eq. 61 becomes 


D 
Sa -r?sin®$+(reosg+5)2 
LL ee 
ise a aay a ‘oe \2C¢ C22 
Ce eS eae 
D 
-r?sin?d+(reosg-5) 2 


+e Cc, * ] 
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Fig. 17. Parallel Nozzles 
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Fig. 18. Radial Profiles of Momentum Flux Density Ratio 
in Two Parallel Jets 
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TABLE 10 
Radial Profiles of Momentum Flux Density Ratio in Two Parallel Jets 
dy = do = 0.898 in.; Ug = 388 ft/sec; D = 3.00 in. 


0.992 0.004 in A 0. 1.652 0.001 1.051 0.006 
0.843 0.036 lls the 0. 1.521 0.004 0.895 0.011 
Or 3s0l= 0.155 OF 1. 0. 1.393 0.006 0.750 0.017 
0.663 0.241 0. 0. 0. L2O2" 0.012 0.635 0.024 
0.599 0.346 0. 0. 0. 0.867 0.037 0.490 0.034 
0.565 0.408 0. 0. 0. 0.736 0.045 0.374 0.043 
0.532 0.432 0. 0. 0. 0.670 0.050 0.248 0.052 
0.500 0.445 Q. 0. 0. 0.601 0.054 OR OOo Td) 
0.470 0.439 0. 0. 0. 0.535 0.056 0.033 0.058 
0.436 0.414 0. 0. 0. 0.475 0.059 
0.425 0.369 0. 0. We 0.410 0.060 0.143 0.057 
0.338 0.266 0. 0. 0. 0.344 0.060 0.263 0.049 
O22 500. 16 0. 0. 0. 0.278 0.060 0.404 0.038 
0.206 0.106 0-263 0.099 0.212 0.059 0.538 0.030 
0.141 0.059 0. 0.197 0.091 0.147 0.059 ONGT2 C2021 
0.075 0.035 0. : 0.131 0.082 0.081 0.058 0.800 0.014 
OF 182202108) +} 0. CROCE 0. 0. 0.917, 0.008 
0.009 0.025 0.248 iY {Ms 6.075 1.061 0.005 
0.057 0.031 0.307 175 Only O09 
0.123 0.055 0. 0. 0. 0.182 0.059 
0.188 0.086 0. 0. 0. 0.248 0.060 
0.254 0.161 Q. 0. 0. 0.370 0.060 
0.310 0.247 0. 0. 0. 0.380 0.061 
0.384 0.333 0. 0. 0. 0.446 0.060 
0.451 0.452 0. 0. 0. Oma e059) 
0.482 0.472 ie 0. 0. 0.574 0.057 
0.517 0.464 kg 0. 0. 0.640 0.054 
0.547 0.432 OB Oo seAg 0.706 0.050 
0.580 0.384 0.565 10, £22 0.969 0.031 
0.646 0.269 0.590 0.118 1.100 0.019 
0.712 0.163 ONT 20 re OF LOZ e232) eOmO el 
0.845 0.046 0.852 0.085 1.363 0.008 
0.975 0.004 0.984 0.055 1.492 0.005 
ig ielale © Oh (Ole: 1.624 0.002 
1.255 0.006 Ale 0. 
1 0. 
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where ¢ is the angle of the radius r measured from the line: 

joining the two nozzles. Rae. 
The momentum flux density ratio along the x-axis 1s 

obtained from Eq. 94 by setting r equal to zero; whence 


a) ees 
Salen py feng ET (95) 


(pw) y/o 


For profiles through the x-axis parallel to the line 
joining the two nozzles, cos. is unity; for profiles 
perpendicular to this line, cos ¢ 1s zero; and for profiles 
along the x-axis, the hyperbolic cosine is unity. These 
equations apply only to the region beyond 10 nozzle diam- 
eters from the nozzles. Furthermore, they involve both the 
assumption of similarity and the principle of superposi- 
tion. For the principle of superposition to be valid, it 
is not necessary that similarity exist, and the principle 
can be tested independently of the similarity assumption by 
the relation 


ik ae: = (pu?) , + (pu7), (96) 


where (pu?), is the experimentally determined momentum 
flux at the point (x,r,¢) when the first jet is operated 
alone, en No is the momentum flux when the second is 
operated alone, and (ou?) , ae 1s the predicted flux when 
both jets are operated simultaneously. 

For comparison, lines representing Eq. 94 with om 
equal to 0.075 are included in Fig. 18. In general, 
the measured values agree with the predicted ones, pro- 
viding support for the principle of superposition. There 
are minor discrepancies. 

Also shown in Fig. 18 is a traverse taken along a radius 
perpendicular to the line joining the nozzles (@ = 90°) ag 
ab. axial’ distance ratio z/d>of 30,.. The solid) lane repre- 
sents Eq. 94, and the agreement is of the same order as was 
obtained by the traverse in the direction of ehe line 
connecting the two nozzles, 

The Reichardt hypothesis together with the principle of 
superposition of momentum flux in the axial direction 
applies well to the flow of parallel jets. Qther measure- 
ments on jets which are not parallel indicate that these 
principles do not apply when the axis of the jets are set 
at an angle so that they impinge on each other. The 
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experimental measurements of x-directed momentum flux differ 
markedly from the values predicted by the superposition of 
momentum flux from the nozzles operating separately as the 
angle between the nozzles is increased. ‘9! 


ce. Single Jet Discharging Parallel to a Plane 


For two parallel jets of equal strength located in the 
Same plane perpendicular to the axis of flow there is a 
plane of symmetry perpendicular to and bisecting the line 
joining the two nozzles. For every increment of x-directed 
momentum that passes through this plane from one nozzle, 
an equal increment passes through in the opposite direction 
from the other nozzle. If an actual plane which is an 
ideal reflector of momentum is placed at this plane of sym- 
metry, then the situation on each side of the plane is the 
Same as before. For a jet discharging parallel to such an 
infinite ideal plane surface the momentum flux distribution 
will be identical with that of one-half the distribution 
from the two parallel jets. 

A physical interpretation of the reflection of momentum 
by solid boundaries is not self-evident. The case was 
investigated experimentally. Equation 94 applies to this 
case, except that the variable D is now to be identified 
with twice the distance between the center of the nozzle 
and the plane. Total-head impact tube measurements were 
made in the flow field of a jet of air discharging at a 
velocity of 388 ft per sec with its axis parallel to a 
large plane surface located 1.50 in. from the axis of the 
nozzle. The inside diameter of the nozzle was 0.898 in. 

The experimental results are shown in Fig. 19 and in 
Table 11, where the momentum flux density ratio is plotted 
against the radial distance ratio r/D. Equation 94 with C, 
equal to 0.075 is represented by the solid lines. The 
agreement is fairly good. The experimental points near the 
plane surface lie above the calculated curves. Beyond the 
maxima the points lie below the calculated curves. A solid 
surface is probably not an ideal reflector of momentum, but 
the analogy is approximately correct. Friction drag 
along the surface does not seem to be important in the 


region near the nozzle. 


19. Applications of Transport in Jets 


Transport in free jets has a number of applications. 
Further measurements are needed to determine the effective- 
ness and the limitations of the methods described. Several 
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Fig. 19 Radial Profiles of Momentum Flux Density Ratio 
in Single Jet Discharged Parallel to Plane Surface 


cases are discussed here as examples of these applications, 
the limitations now known to exist are indicated, and 
the measurements needed as a basis for confirming the 
method are mentioned, 

Consider a case where smal] droplets of water are formed 
in a stream of hot air and these spread out into the stream 
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Table |} 
Radial Profiles of Momentum Flux Density Ratio 
in Single Jet Discharged Parallel to Plane Surface 
d = 0.898 in.; Up = 388 ft/sec; D = 3.00 in. 


0. 0.045 0. 0. 0. 0. 0.057 0.049 
0. 0.050 0. 0. Q. 0. O59 OAs) 
0. 0.056 0. 0. 0. 0. 0.129 0.064 
0. 0.068 0. 0. 0. 0. 0.162 0.066 
0. 0. 089 0. Or 0. 0. 0.195 0.068 
Os, Oedialy 0. 0. 0. OF 0.228 0.068 
0. 0.145 0. 0. 0. 0. 0.261 0.067 
0. Oss) 0. 0. 0. 0. 0.359 0.066 
0. 0.224 0. 0. 0. 0. 0.425 0.063 
0. OE 269 0. 0. 0. 0. 0.470 0.059 
0. 0.332 0. 0. OF 0. 0.491 0.058 
OF 0.362 0. 0. 0. 0. 0.622" 0.049 
0. 0.425 0. 0. 0. 0. Osis WSO 
0. 0.458 OF 0. OF 0. 0.885 0.029 
0. 0. 467 O}; 0. 0. ‘OF LOL 0017 
0. 0.419 IL. 0. ile 0. 1.143 0.009 
OF 0.301 i; Or Ale, 0. 2 303006 
Oe 0. 083 iE 0. 

0. 0. 

Li: 0. 


as a free jet. Heat is transferred from the hot air to the 
surface of the droplets; at the same time water vapor 
leaves the drops and is transported into the stream. When 
the drops are very small they follow the fluctuations of 
the air stream. It is then possible to write the differen- 
tial equations for the transport of water vapor, liquid 
droplets, heat, and x-directed momentum flux. Water vapor 
is added to the air stream by evaporation while heat is 
removed, and terms are required in the equations to account 
for this addition and removal. These differential equations 
may be modified according to the generalized Reichardt 
hypothesis and solved by numerical methods. The solution 
gives the distribution of temperature and water vapor 
concentration and also the rate of evaporation of the water 
droplets as they move downstream in, the rapidly changing 
environment of the stream. Measurement of the distribution 
of temperature and water vapor concentration at several 
sections downstream from the nozzle would serve for compar- 
ison with the calculated distribution. 
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Another application is the free jet turbulent diffusion) 
flame. Here a jet of fuel gas issues into air, entraining} 
and mixing with the air. The gas and air intermingle, 
diffusing into one another. At the proper concentration| 
and temperature a chemical reaction occurs. The pro-: 
ducts of the reaction and the heat from the reaction 
are transported throughout the flame. This case is compli-: 
cated by the rate of the chemical reaction or by the: 
discontinuities in concentration and temperature when, 
the reaction is rapid. 

The coalescence of small droplets is important in 
atomizing and in the handling of aerosols. Methods such as: 
those described above may be used to calculate the trans- 
port of aerosol particles issuing from a small source and 
spreading out in a free jet. By tracing the concentration 
and rate of coalescence throughout the jet, the overall 
effect of coalescence may be calculated. 

A jet of gas issuing into a tube or duct and inducing 
the flow of gas through the duct by entrainment is related 
to the free jet. The momentum flux distribution across the 
duct at several sections was measured as well as the static 
pressure at the wall along the length of the duct. It was 
found that the static pressure changes rapidly in the 
vicinity of the jet source; in this respect the situation 
differs from a free jet where the static pressure is 
essentially uniform. However the momentum when combined 
with the pressure force in a single term permits extension 
of these methods to the jet in a duct. 

Interesting qualitative conclusions have resulted from 
the measurements in a duct. The changes in static pressure 
along the length of the duct in the vicinity of the jet re- 
sult principally from the redistribution of x-momentum flux 
density over successive sections of the Stream. The pres- 
Sure at first rises to a maximum and then decreases. These 
pressure changes differ from those observed in fully 
developed flow through a pipe. The momentum tends to become 
uniformly distributed over a section, so that a relatively 
flat velocity profile may be observed a few duct diameters 
downstream from the jet: sthis profile then gradually 
becomes more rounded and approaches the normal distribution 
in fully developed turbulent flow. 

The action of a free jet in entraining the surrounding 
fluid-is <a natural function, and if fluid to meet this 
entrainment requirement is not supplied the jet becomes un- 
stable. This instability was readily observed with the jet 
in the duct when the supply of secondary air for entrain- 
ment was progressively reduced. With ample Secondary air 
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the jet was stable and the flow chiefly in one direction. 
As the supply of secondary air decreased, the entrainment 
air was supplied from further downstream by sucking the air 
“hext to the walls back to the jet. Thus a reversed flow 
peddy is set up, and this is a condition of instability. 
The same pattern is evident in the region downstream from 
an orifice plate in a pipe where the supply of fluid for 
entrainment in the jet issuing from the orifice is entirely 
restricted. Similar eddies result from the entrainment 
demand when a stream flows past a baffle. 


V. SUMMARY 


A method based on Reichardt’s hypothesis of momentum 
transfer is described, and is shown to be useful ang 
accurate for complex cases of the transport of heat, 
mass, and momentum in free jets within the limits of 
precision of the experimental data. The method is not 
generally applicable to ail cases of turbulent transfer, 
but it does permit the solution of many special problems 
which were hitherto unsolved. 

The theory is developed from the basic equations of 
motion and continuity. The physical models proposed by 
Prandt] and extended by Taylor were also examined; however 
the hypothesis proposed by Reichardt was selected as offer- 
ing greater promise of useful application because it 
changes the nonlinear flow equations to linear differential 
equations and simplifies the solution of complex cases. 
Peichardt’s hypothesis states (Eq. 49) that the mean of the 
product of the instantaneous axial and radial velocities 
1S proportional to the change in the mean of the square of 
the axial velocity with distance along a radius. : 

The extension of Peichardt’s hypothesis to mass and heat 
transport is proposed. This more general statement, 
contained in Eq. 53, is that the mean of the product of the 
radial velocity and a transport quantity is proportional to 
the change of the mean of the product of the axial velocity 
and the transport quantity with distance along a radius. 
The statement is substantiated by experiment when the 
radial’ velocity ‘is small relative touthesasial velocity, 
but is not valid for large radial velocities. 

The extended hypothesis is compared with experimental 
data from single free jets covering the transport of momen- 
tum, heat and mass,’ and of: momentum, ini the complex region 
near a nozzle. The transport of momentum from multiple 
parallel jets and from a jet issuing parallel to a wall is 
described. The applications to evaporation from a spray of 
liquid droplets in a jet and to combustion in a similar 
case are discussed. These principles also provide a basis 
for examining the coalescence of particles in a turbulent 
jet and hold promise of being applicable to the spreading 
of a jet issuing into a duct when the flow of the secondary 
stream in the duct is sufficient to stabilize the Neu 
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